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Detection of ultra-high energy neutrinos will be useful for unraveling the dynamics of the most vio¬ 
lent sources in the cosmos and for revealing the neutrino cross-section at extreme energy. If there ex¬ 
ists a Greisen-Zatsepin-Kuz’min (GZK) suppression of cosmic-ray events above Eqzk ^ 5 x 10*^® eV, 
as predicted by theory, then the only messengers of energies beyond Agzk are neutrinos. Gosmic 
neutrino fluxes can initiate air-showers through interaction in the atmosphere, or in the Earth. Neu¬ 
trino trajectories will be downgoing to nearly horizontal in the former case, and “Earth-skimming” 
in the latter case. Thus it is important to know the acceptances (event rate/flux) of proposed 
air-shower experiments for detecting both types of neutrino-initiated events. We calculate these 
acceptances for fluorescence detectors, both space-based as with the EUSO and OWL proposals, 
and ground-based, as with Auger, HiRes and Telescope Array. The neutrino cross-section is 
unknown at energies above 5.2 x 10^® eV. Although the popular QCD extrapolation of lower-energy 
physics offers the cross-section value of 0.54 x 10“®^ (i5^/10®®eV)® ®® cm®, new physics could raise or 
lower this value. Therefore, we present the acceptances of horizontal (HAS) and upgoing (UAS) air 
showers as a function of over the range 10“®“^ to 10“®® cm®. The dependences of acceptances 
on neutrino energy, shower-threshold energy, shower length, and shower column density are also 
studied. We introduce a cloud layer, and study its effect on rates as viewed from space and from 
the ground. Eor UAS, we present acceptances for events over land (rock), and over the ocean (wa¬ 
ter). Acceptances over water are larger by about an order of magnitude, thus favoring space-based 
detectors. We revisit the idea of Ref. Q to infer at E^, ^ 10®® from the ratio of HAS-to-UAS 
events, and obtain favorable results. Included in our UAS calculations are realistic energy-losses 
for taus, and Earth-curvature effects. Most of our calculation is analytic, allowing insight into the 
various subprocesses that collectively turn an incident neutrino into an observable shower. 

PACS numbers: 13.85.Tp;95.55.Vj;95.85.Ry;96.40.Pq 


I. INTRODUCTION AND PURPOSE 

Detection of ultra-high energy [E^, > 10^® eV = EeV) neutrinos is important for several reasons. First of all, 
neutrino primaries are not deflected by magnetic fields and so should point back to their cosmic sources. This contrasts 
with cosmic-rays, which are charged and follow bent trajectories. Secondly, well above the Greisen-Zatsepin-Kuzmin 
(GZK) energy of Eqzk ~ 5 x 10^® eV is, they may be the only propagating primaries. As such, they may be the 
only messengers revealing the ultimate energy-reach of extreme cosmic accelerators, generally believed to be powered 
by black holes. Above Eqzk, the GZK suppression dlllii of cosmic-rays results from the resonant process 
^ + 7CMB ^ A ^ A -b tt; Eqzk is the lab-frame energy corresponding to the kinematic threshold -/s = Ma for 
excitation of the intermediate A resonance. A handful of cosmic-ray events have been detected with estimated energies 
exceeding I0^° eV. The record energy is the famous Fly’s Eye event at 3x I0^° eV The observable neutrino spectrum 
could extend to much higher energies. Thirdly, in contrast to cosmic-rays and photons, neutrinos are little affected 
by the ambient matter surrounding the central engines of Nature’s extreme accelerators. Accordingly, neutrinos may 
carry information about the central engine itself, inaccessible with other primaries. In principle, neutrinos may be 
emitted from close to the black hole horizon, subject only to energy-loss due to gravitational redshifting. An analogy 
can be made to solar studies performed with photons versus neutrinos. The photons are emitted from the outer 
centimeter of the Sun’s chromosphere, while the neutrinos are emitted from the central core where fusion powers the 
Sun. Fourthly, neutrinos carry a quantum number that cosmic-rays and photons do not have - flavor. Neutrinos 
come in electron, muon, and tau flavors. One may think of this “extra” flavor degree of information as the neutrino’s 
superb analog to polarization for the photon, or nucleon number A for the cosmic-ray. Each of these attributes, flavor. 
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polarization, and nucleon number, carries information about the nature and dynamics of the source, and about the 
environment and pathlength of the inter-galactic journey. The flavor ratios of cosmic neutrinos are observable Q]- 
Several papers have recently analyzed the benefits that neutrino flavor identification offers for unraveling the dynamics 
of cosmic sources |^. The fifth reason why ultra-high energy neutrino primaries traveling over cosmic distances are 
interesting is that such travel allows studies of the fundamental properties of neutrinos themselves. For studying some 
properties of the neutrino, such as neutrino stability/lifetime pseudo-Dirac mass patterns M , it is the cosmic 
distance that is essential; for other properties, it is the extreme energy that is essential. A clear example of the latter 
is any attempt to determine the neutrino cross-section at energies beyond the reach of our terrestrial accelerators. 

In this paper we will examine the potential for cosmic-ray experiments designed to track ultra-high energy air- 
showers by monitoring their fluorescence yield 0. to detect horizontal air-showers (HAS) and upgoing air-showers 
(UAS) induced by a cosmic neutrino flux. We will also study the ability of these experiments to infer the neutrino- 
nucleon cross-section a^N at energies above 10^® eV, from the ratio of their UAS and HAS events. Such energies are 
orders of magnitude beyond the energies accessible to man-made terrestrial accelerators. From the point of view of 
QCD, such a cross-section measurement would be an interesting microscope into the world of small-x parton evolution. 
The neutrino cross-section above 10^® eV could agree with any of the various QCD-motivated extrapolations that 
have been published nnn, or not. The cross-section could also be quite different than the extrapolations. For 
example, if a new threshold is crossed between terrestrial neutrino energies ~ 100 GeV, and the extreme energies 
reached by cosmic-rays, ^ 10^^ GeV, then the cross-section could much exceed the QGD-extrapolations. On the 
other hand, saturation effects can significantly reduce the total cross-section at these very high energies ^2. The 
nine orders of magnitude increase in lab energy reach corresponds to 4.5 orders of magnitude increase in center-of- 
momentum energy reach. Even the center-of-momentum energy at the e — p HERA collider is more than three orders 
of magnitude below the cosmic-ray reach. This remarkable energy reach of cosmic-rays presents ample room for new 
physics beyond our Standard Model. Proposals for new physics thresholds in this energy region include low-scale 
unification with gravity, in which neutrino-nucleon scattering produces mini-black holes |l5l| and/or brane-wraps m, 
non-perturbative electroweak instanton effects 0 , compositeness models 0 , a low energy unification scale in string 
inspired models Kaluza-Klein modes from compactified extra dimensions . All of these models produce a 

strongly-interacting neutrino cross-section above the new threshold. Dispersion relations allow one to use low-energy 
elastic scattering to place constraints on the high-energy cross-section [^, but the constraints are quite weak. 

For HAS and UAS, we provide analytical calculations of the event-rate to flux ratio as a function of a^N- This 
ratio is known as the “instantaneous experimental acceptance”, with units of areaxsolid angle. The time-averaged 
acceptance includes an experimental “duty factor,” the fraction of time that the experiment is functioning. We will 
not include the duty factor in our calculations of acceptances. We note that acceptances are also sometimes called 
“apertures.” 

Experimental acceptance offers a very meaningful figure of merit for statistical reach. One has merely to multiply 
an experiment’s acceptance by Nature’s flux to arrive at an event rate for the experiment. Multiplying again by the 
experiment’s run time (including the duty factor), one obtains the total number of events. Acceptance times run time 
is termed the experimental “exposure”. 

The acceptances we calculate are scalable to large area experiments such as HiRes, Auger, and in the near future 
Telescope Array, which are anchored to the ground, and to super-large area experiments such as EUSO and OWL, 
which are proposed to orbit the Earth from space. A horizontal shower, deeply initiated, is the classic signature for 
a neutrino primary. The weak nature of the neutrino cross-section means that horizontal events begin where the 
atmospheric target is most dense, low in the atmosphere. In contrast, the ultra-high energy pp cross-section exceeds 
100 mb, so the air-nucleon cross-section exceeds a barn! Even the vertical atmospheric column density provides 
hundreds of interaction lengths for a nucleon, and so the cosmic-ray interacts high in the atmosphere. The weak 
nature of the neutrino cross-section also means that the event rate for neutrino-induced HAS is proportional to the 
neutrino-nucleon cross-section. 

For an neutrino-induced UAS, the dependence on neutrino cross-section is more complicated, and more interesting. 
The Earth itself is opaque for neutrinos with energies exceeding about a PeV of energy. However, “Earth-skimming” 
neutrinos, those with a short enough chord length through the Earth, will penetrate and exit, or penetrate and 
interact. In particular, there is much interest in the Earth-skimming process > r in the shallow Earth, followed 
by T decay in the atmosphere to produce an observable shower. In Ref. Q it was shown that the rate for the Earth- 
skimming process i>r t is inversely proportional to cr^Ar. There it was emphasized that could be inferred from 
a measurement of the ratio of HAS to UAS rates. ^ Of course, an implicit assumption is that there is enough neutrino 


^ The prospects of inferring the neu trino-nucleon cross-section in the energy range of 100 TeV - 100 PeV at neutrino telescopes such as 
IceCube, were studied in Ref. [2^; prospects at higher energies were studied in Ref. |2^ for the Auger observatory. 
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flux at extreme energies to generate HAS and UAS event samples. 

The inverse dependence of UAS rate on is broken by the r ^ shower process in the atmosphere. As the 
cross-section decreases, the allowed chord length in the Earth increases, and the tau emerges with a larger angle from 
the Earth’s tangent plane. This in turn provides a smaller path-length in air in which the tau may decay and the 
resulting shower may evolve. This effect somewhat mitigates the inverse dependence of the UAS on 

Ref. m provided an approximate calculation of the whole UAS process, and gave an approximate result for the 
dependence of the HAS/UAS ratio on cTyAr. In this work, we improve upon Ref. Q in several ways. We include the 
energy dependences of the tau energy-losses in the Earth, and of the tau lifetime in the atmosphere. For the energy- 
losses, we distinguish between tau propagation in earth rock and propagation in ocean water. These calculations are 
carried out in Section QH On the issue of shower development, we incorporate the dependence of atmospheric density 
on altitude. We also impose requirements on the resulting shower such that a sufficiently long visible shower-length 
is projected onto the Earth’s tangent plane, thus meeting experimental requirements for visibility. This is done in 
Section ITTH In the case of the upgoing showers, the pathlength of the pre-decayed tau may be so long that the Earth’s 
curvature enters into the altitude dependence. We include the non-negligible correction from curvature in Section ITm 
We include the partial loss of visibility due to high cirrus or low cumulus cloud layers in Section HU It is estimated 
that clouds will obscure the viewing area about 60-70% of the time. For ground-based observation, it is mainly the 
low-lying cumulus clouds that limit visibility. For space-based observation, it is mainly the high cirrus clouds that 
limit visibility. ^ And in Section El we combine the corrections from clouds with that from the Earth’s curvature. 

Our results are illustrated in a series of plots of acceptances, for ground-based and space-based experiments, versus 
neutrino-nucleon cross-section, in Section ED Situations with and without cloud layers are analyzed, as are events 
over solid earth and over the ocean. Incident neutrino energies, energy thresholds for experimental detection of 
the air-shower, and various shower-trigger parameters are varied. Earth-curvature effects are included in our UAS 
calculations. These reduce the event rate. Next comes the discussion Section EID It presents several small issues, 
and includes a comparison of our work with prior work. A final Section recaps our conclusions. Some of the more 
tedious but necessary formulas are derived in an Appendix. 

The reader who believes that a picture (or four) is a worth thousand words may wish to jump to Section l\TI Such 
a reader especially may find it useful to reference Tables ^ and m where the variables and parameters are defined. 

Among our conclusions, we find that the HAS/UAS ratio is or order of unity for cross-section values very near to 
the commonly extrapolated value of 0.5 x 10“^^ cm^ at Ei, ~ I0^° eV. This is fortunate, for it offers the best possibility 
that both HAS and UAS rates can be measured, and a true cross-section inferred. We display our HAS and UAS 
acceptance plots for a cross-section range from superweak 10“^"^ cm^ to a microbarn, cm^. This range includes 

the QCD-extrapolations of and the region of the HAS/UAS cross-over. It also encompasses any effects of new 
neutrino physics, either increasing or decreasing cr^^. The highest energy for which the neutrino cross-section has 
been measured is that at the HERA accelerator. The measurement is ^ 2 x 10“^‘*cm^ at yG = 314 GeV 1^ . 
the latter corresponding to an energy on fixed target of 5.2 x 10^^ eV (52 TeV). It is hard to imagine that at 
10^° eV would not have grown beyond the HERA value. Even so, the acceptances shown for superweak cross-sections 
may have some relevance to a possible WIMP flux . Modeling of a WIMP event rate requires modifications in the 
shower development for HAS, and in the chain WIMP^UAS, that we do not pursue here. 


II. AIR-SHOWER RATES 

The variables and parameters needed to describe UAS and HAS are sufficiently numerous that we have collected 
many of them in Table HI and ITU for easy reference. In Table ITTII we explain the different symbols used throughout this 
work. Many of the variables are best explained by the three schematic diagrams in Figs. [21-El 


A. Upgoing Air-Showers (UAS) 

Ultra-high energy (UHE) neutrinos are expected to arise from the decay of pions and subsequently muons produced 
in astrophysical sources Q. For this decay chain, the flavor mix at the source is Ve ■ >^ii = ^ ‘2- The maximal mixing 
between and Vr inferred from terrestrial oscillation experiments then leads, after propagation for many oscillations 


^ In fact, low-lying cumulus clouds may aid in HAS identification for space-based observing. When the HAS hits the cloud layer, diffuse 
reflection of the forward Cerenkov cone can be seen as a one-time “Cerenkov flash”. The time of the flash and the measured height of 
the cloud then provide the absolute (t,z) coordinates of the shower. 
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L 

chord length of v trajectory through Earth 

Zint 

vertical height (depth) of HAS (UAS) v interaction 

-2^dk 

altitude of upgoing Vt decay (no Earth-curvature) 

^dk 

altitude of upgoing v-r decay including Earth-curvature 

zu 

maximum visible shower altitude (HAS ^ UAS) 

Zl 

minimum visible shower altitude (HAS yf UAS) 

Zb 

altitude where shower first attains threshold brightness (HAS yf UAS) 

Ze 

altitude where shower extinguishes (HAS yf UAS) 

crit 

“^cloud 

critical altitude for suppression from cloud layer 

4(UAS) 

altitude where UAS attains threshold brightness, including Earth-curvature 

4(UAS) 

altitude where UAS extinguishes, including Earth-curvature 

e. 

zenith angle of HAS event 

dn 

nadir angle of UAS event (no Earth-curvature) 

^*hor = ^ — 02 

horizontal angle of UAS event 

e'r. 

nadir angle of UAS event including Earth-curvature 

^hor 

horizontal angle of UAS event including Earth-curvature 

dtot 

total column along chord of Earth 

di/ 

column density of v in the Earth 

d-T 

column density of r in the Earth 

cos 6s 

minimum shower angle, cloud-dependent, for space-observatory 

cos 9q 

minimum shower angle, cloud-dependent, for ground-observatory 

2(HAS) 

maximum altitude from which initiated HAS can reach the ground 

2(UAS) 

minimum altitude from which initiated UAS can reach Zthin 

2'(UAS) 

5(UAS) with Earth-curvature included 

2(UAS (g) G) 

5(UAS) modified for cloud layer above 

2(HAS (g) S) 

£(HAS) modified for cloud layer below 


TABLE I: List of variables and their meaning. (Conversion between variables « and w is given by wcosO = z.) 


h 

scale height of the atmosphere 

8 km 

Aground 

ground altitude, kept as a symbol for later substitutions 

0 

-2^thin 

altitude beyond which air is too thin to fluoresce significantly 

3h 

/3i9 

tau energy-attenuation constant at Et = 10^® eV 

1.0 (0.55) xl0“®cm^/g for rock (water) 

a 

exponent of the energy-dependence of /3 t 

0.2 

d-vert 

vertical atmospheric column density 

1,030 g/cm^ 

C^hor 

horizontal atmospheric column density 

36,100 g/cm^ 

C^min 

minimum acceptable shower column density 

300, 400 g/cm^ 

duxax. 

maximum shower column density at extinction 

1200 , 1500 g/cm^ 

^min 

minimum acceptable shower length projected on the Earth’s surface 

10 km, 5 km 

Rfov 

radius (or half-scale) of the experimental field of view 

230 km 

2w 

depth of ocean 

3.5 km 

'2'cloud 

altitude of cloud layer 

2 , 4, 8, 12 km 

E, 

incident neutrino energy 

10 ^°, 10^^ eV; 

Tpsh 

-^th 

detector threshold energy 

10^®, 5 X 10^® 

sCh 

tau threshold energy 

f (3) X for hadron (electron) mode 


neutrino (or WIMP) cross-section 

10"®°, 10 ^®\ 10"“, 10"®® cm® 


TABLE II: List of parameters, their meaning, and their chosen value(s); the bold-faced value is the chosen “canonical” value. 





5 


HAS 

Horizontal Air-Shower 

HAS 

Upgoing (“Earth-skimming”) Air-Shower 

HAS0S 

HAS seen from space-based observatory 

UAS®S 

UAS seen from space-based observatory 

HAS(g)G 

HAS seen from ground-based observatory 

UAS(g)G 

UAS seen from ground-based observatory 


TABLE III: List of symbols and their meaning. 


lengths and to a very good approximation, to a flavor ratio at Earth oi Vg : : Vr = ^ ^ i-e., flavor democracy. 

Thus, a detector optimized for iZg or or v-r can expect a healthy signal from cosmic neutrinos. 

It is useful to define the neutrino charged-current (CC) interaction mean free path (MFP) as 


A. 


1 63 km 

PP _ — 

P ^'2.65 0-31 


( 1 ) 


where p is the mean number-density of the target-matter, and p 2.65 is the mean density in units of the value for 
surface rock psr = 2.65 g/cuY. Density is usually expressed in units of g/cm^, with the multiplicative factor of 
IVa = 6.022 X 10^^ g“^ implicitly understood. The mean density of ocean water is 1.0 g/cm^. The cross-section 0-31 
is the CC cross-section in units of 10“^^ cm^. The commonly used high-energy neutrino-nucleon CC cross-section 
extrapolated from QCD is 0.54 x 10“^^ (ify/10^°eV)°’^®^cm^. 

We will ignore the NC contribution to the neutrino MFP for three reasons. First, the NC cross-section is expected 
to be small compared to the CC cross-section, as it is known to be at the lower energies of terrestrial accelerators. 
Secondly, the NC interaction does not absorb the neutrino, but rather lowers the energy of the propagating neutrino 
by a small amount; in the SM, the energy loss is only {y) ^ 20%. And thirdly, the increase in complexity of our 
calculation, when the NC MFP is included, seems unwarranted. We also ignore multiple CC interactions due to the 
“tau regeneration” decay chain ^ t ^ v^. Here, it is the long decay length of the tau that results from production 
at Eh > 10^^ eV that makes tau regeneration negligible. 

In Fig. n we show an interesting relation between the neutrino cross-section, the neutrino’s MFP in the Earth, and 
roughly speaking, the maximum horizontal angle for which the neutrino may transit the Earth. 

In this figure, the Earth has been approximated according to the two-shell model. There is a central core with mean 
density 12 g/cm^ out to a radius of 3485.7 km, and a mantle with mean density 4.0 g/cvo? out to the the Earth’s 
radius of = 6371 km. The point of this figure is that, although the Earth is marginally transparent for neutrinos 
with the HERA cross-section of 2 x 10“^^ cm^, the Earth quickly becomes opaque at larger cross-section. For the 
cross-section values extrapolated to ~ 10^° eV, horizontal angles are very small, and the trajectories are truly “Earth 
skimming” (cf. Eq. 0). 

The tau energy-attenuation length is At- = (/?,-p)“^, where Pt{E) is the coefficient giving a scale to tau energy loss: 


dEr/dx = -Pr{E)pEr. (2) 

The coefficient f3r{E) is weakly energy-dependent. For the energies of interest, E^ > 10^® eV, tau energy losses 
are dominated by photo-nuclear pro cesses, being the electromagnetic mechanisms of ionization, bremsstrahlung and 
e^-e” pair production negligible |26l|. We find that the recent calculations of l3r{E) 113, are well fitted in the 

energy region of interest by a simple power law |^: ® 

/3-(^)=/3i9 > « = 0.2, (3) 


with the constant pre-factor /3®g = l.Ox 10“® cm^/g for surface rock ((H) = 22, {Z) = 11), and PYg = 0.55x 10“® cm^/g 
for water ((H) = 11.9, {Z) = 6 . 6 ); Pt{E) scales as (H). The tau energy-attenuation length at Er = 10^® eV is 
At = 3.8 km in surface rock, and 18 km in water. The tau decay MFP is ctt = 490 (Et/IO^^cV) km. For taus with 
energies at and above 10^® eV, the decay MFPs are much longer than the energy-attenuation length. In this paper. 


^ In a very recent paper m, a logarithmic fit to Pt{E) is presented. We find that our fit agrees quite well with that one in the region of 
our interest, 10^® eV < Et/ < 10^^ eV. For our purposes, the power law fit is more useful in that it allows analytic integration of some 
energy-dependences. 
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FIG. 1: Shown are neutrino trajectories for which the interaction MFP matches the chord length through the Earth. The 
various trajectories are parameterized values of the neutrino cross-section. Also shown is the trajectory’s angle with respect to 
horizontal. 


we safely neglect the small probability of decay within the Earth for those tans which would otherwise emerge from 
the Earth with energy above 10^® eV. We have checked that the results we present in this work are reduced by less 
than a few percent when the tau decay probability within the Earth is included. 

The muon energy-attenuation length is 7 times smaller than that of the tau, and the electron energy-attenuation 
length is many times smaller again (the /x decay length is ~ 10® times longer than that of the tau). Because the 
energy-attenuation length for a tau is an order of magnitude longer than that of a muon, UAS events are dominantly 
initiated by the CC interaction of tau neutrinos. 

The ratio of the tau energy-attenuation length to the neutrino MFP Xt/^v = NAcr^^/f^r ~ (cr^^/10“®^cm^) x 
0.06 (0.11) for rock (water), is independent of p and only weakly dependent on tau energy. For ^ 2 x 10“®®cm^, 
we expect most of the path length in Earth (rock or water) to be neutrino; for ^ 2 x 10“®°cm^, we expect most 
of the path length to be tau. In detail, this remark will also depend on the direction of the initial neutrino through 
the total chord length, and on the threshold energy of the detector (the minimum record-able tau energy). 

Consider a tau produced in the Earth along a chord on the trajectory of an incoming Label the chord length 
by L and the distance between the interaction and the Earth’s surface by Wint, as shown in Fig. [21 In general, we 
will use the variable w to represent distance or location along the lepton trajectory; when needed, z will label the 
distance normal to the Earth’s surface, and x and y will label the Earth’s tangent plane. The chord length L and the 
nadir-angle of the upgoing neutrino trajectory characterize the same degree of freedom, coupled together by the 
geometric relation L = 2i?0 cos0„, where i?® = 6371 km is the radius of the Earth. 

We follow the calculation of the rate of UAS events as given in Ref. Q, although with some important improvements. 
The rate is 

TAT 

(UAS) = F.^ttA / ^ P.^^AL) 0 Pdk(L), (4) 

Jo 

where Pu^^t{L) is the probability for a v^. along a chord of length L in Earth to produce a tau which exits the surface 
with a given threshold energy Uth, Tdk(L) is the probability of decay for the tau emerging into the air, and is 
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z 



the Vr differential flux in the usual dimensional units of (energy • time • area • steradian) The ® symbol denotes 
coupling between the two probabilities, as shown in detail below. In Eq. 0 we have not considered the possibility 
of the tau decaying inside the Earth, which for the energies of interest in this study is a very good approximation. 
The detector’s field of view (FOV) is >1 = J dx J dy. In general, the FOV includes some detector-efficiency weighting 
for shower identification, as explained in Section uni The operator ttA f LdL/2R^ is a convenient rewriting of the 
integrals / dA ■ h f dfl, obtained when use is made of the relation L = 2 i ?0 cos0„. The angular dependence of the 
interaction and decay probabilities are therefore implicit in the L-dependence. 

The interaction probability is given by 




^min{L ,iuth} 


dw 






(5) 


where dwint/^v is the probability for neutrino conversion into tau lepton in the interval [rcint,iCint + dti'int], and the 
exponential gives the survival probability of the neutrino to reach the interaction point Wint- The column density 
traversed by the neutrino is given by 


d^= dw/dearth(w) ■ (6) 

'ii’int 

If the density Pearth were a constant, the exponential in Eq. © would be simply e ^ i _ ^CC 

Such will be the case if the absorption in the Earth limits chord lengths to just the outer layer of Earth-matter (water 
or surface rock). The angle of the trajectory above the horizon is related to the chord length as sindhor = L/2Rq. 
Setting the chord length equal to the neutrino MFP A^, we get for the typical angle 

dhor (2 i?© Pearth)-^ = 0.28° X agY . (7) 

A commonly quoted extrapolation for the neutrino CC cross-section is (T 31 = 0.55 at 10^° eV |l3- Comparisons 
with the “critical” angles delimiting the various density boundaries in the Earth, given in Table © in the Appendix, 
then reveals that over ocean, ^ 4 x 10“^^cm^ gives rise to events whose trajectories were dominantly in only 
water; and over land, ^ 10 “^^cm^ gives rise to events whose trajectories were dominantly in only surface rock (as 
opposed to mantle or core). For these events, the Earth density is approximately a constant. We also study smaller 
cross-sections, for which the density is not constant along the path-integral. In the Appendix we present our general 
calculation of d^. 

The bound Wth(drt) on the depth of Wint integration in Eq. lOl is determined by the requirement that the tau emerge 
from the Earth with sufficient energy, to produce air showers which trigger the detector apparatus. In general, 
Wth is angle-dependent because the density in the Earth is angle-dependent. The mean energy of the tau emerging 




( 8 ) 


from the Earth is obtained by integrating Eq. m. The result is 

Eo 


Er (u?int) — 


l + aX(ri;int) 


-| 1 /C 


where Eq = (1 — (y)) is the mean energy of a tau created by an incoming neutrino with incident energy E^, and 
{y) is the average inelasticity parameter which we will take as {y) = 0.2 Thus, we take Eq = Q.SE^. We 

define T{wint) as the dimensionless tau “opacity” from point of production to the Earth’s surface, normalized to a tau 
with E = IQi® eV, 




dw I3 i9{z) Pea.ith{z) ■ 


(9) 


This definition allows isolation of the energy-dependence of Pt{E) in a separate factor, evident in Eq. ®. Note that 
both /3 i9 and pearth in m depend on the Earth’s composition (e.g., water versus rock), which in general changes 
with depth z. For UAS rising from land, there is no z-dependence: Pearth = Psr and Pig = 1.0 x 10“®cm^/g are 
fixed, and the tau opacity is simply the path length in units of the tau energy-attenuation length at Eq = 10^® eV, 
^Cint/Ar(10^®eV), with Ar(10^®eV) = (/JigPsr)”^- On the other hand, for UAS rising from the oceans, there is a 
discontinuity at the ocean’s bottom: Pearth comprises two contributions, one from ocean water and the other from the 
underlying rock. In the Appendix, we show the calculation of X for this case. We will take the depth of the ocean Zw 
to be 3.5 km. 

Setting Eq. © equal to E^^^, one obtains 


X{wth{dn)) = - 
a 


lO^^eV 


th 


/l0i®eVy 

J 


( 10 ) 


as the equation defining the integration maximum icth- ^ For UAS rising from land, Pig psr is constant and the 
integration and inversion of Eq. to get wth is trivial. For UAS rising from the ocean, the integration and 
inversion of Eq. m to get Wth is more complicated, as the path comprises a water and a rock component. Both 
cases, land and ocean, are dealt with in the Appendix. 

The decay probability is 


F’dk(F) = 


f 


-e 


with the tau lifetime in the lab-frame given by 


Er 


■ Trf = 


392 (E^/lO^^eV) km 


[l + aX(wint)(0.8F;,/10i9eV) 


a I a 


( 11 ) 


( 12 ) 


where trf is the rest-frame value of the tau lifetime. The numerical expression in Eq. m properly includes the 0.8 
mean factor for energy transfer between the incident Vr and the r. We remind the reader that for the simple case of 
UAS over rock, the opacity is just X = Pig psr Wint. The more complicated case for UAS over oceans is dealt with in 
the Appendix. 

When the tau decays, it has a 64% branching probability to decay to Vr + hadrons. For an unpolarized tau, ~ 2/3 
of the energy goes into the hadrons, and therefore into the shower. Accordingly, for this mode we take the relation 
between tau and shower energies to be We define A®)) to be the minimum-energy trigger for the detector. 

Thus, we have the threshold relation E^^ = | E/J). The tau also has 18% branching probabilities each into n + v + e 
and n + D + p. The electronic mode immediately creates an electromagnetic shower with ^ 1/3 of the tau energy, on 
average. So for the electronic mode we take the relation between thresholds to be = 3E®)(. The muonic mode 
is ignorable, for the decay length of the muon exceeds the distance to the ground. In our calculation of the UAS 
acceptance, we will weight each tau decay with 64% for the hadron mode where i?®)) = | and 18% for the electron 
mode where U®)) = A E^^^; the remaining 18% is the unobservable muon mode. ® 


^ The a —S' 0 limit of the In of the denominator in Eq. is easily seen to be X, and so the q ^ 0 limit of Er is Eq e ^. 

^ Writing and similarly for one readily finds the q —> 0 limit of Eq. 1101 to be X{wth{0z)) = In (Eq/E^^). 

® The tau is 100% polarized (to order mr jEr') at production in the Earth. It is possible that even after multi-scattering in the Earth 
(mainly due to photo-nuclear interactions), the tau retains some of its initial polarization. If so, then the decay particle with helicity 
opposite to that of the tau is softer on average. The net result is slightly more energy transferred to the electromagnetic shower, and 
slightly less energy transferred to the hadronic shower in. 
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FIG. 3: Lateral snapshot of UAS and HAS; 2 : labels vertical altitudes and depths. 


The two integrals in Eqs. © and are coupled via the tCint-dependent lifetime of the tau. When we later introduce 
constraints due to cloud covering, we will see further coupling among the integration variables. The exponential in 
Eq. (Hill describes the survival probability of the upgoing tau lepton to reach the decay point Wdk- There is some 
probability for the tau to decay inside the Earth in the interval [0,Wint], but as we mentioned above it is negligibly 
small in the energy range of interest. There is regeneration of tau neutrinos over the whole Earth due to the tau 
production and decay chain, but the regenerated taus with their lower energy contribute negligibly to the high-energy 
sample discussed here and so are not included. 

In practice, a sufhcient column density of air beyond the tau decay point Wdk is required such that the decay products 
fully develop into a shower. This requirement will cutoff the integration in Eq. m, and provide an L-dependence 
(or cos 0„-dependence) to Pdk- In the original study 0, a simple analytic result for the decay integral was 
obtained by invoking certain approximations. The tau lifetime was taken to be a constant over the energy-range of 
interest, and the integral was cutoff at the scale height of the atmosphere, h = 8 km. With these approximations, one 
obtains for the decay integral Pdk = 1 — = 1 — Also, the air shower rate per incident Vr was 

computed analytically for the case where the angle above the horizon satisfies ddor ^ (lO^^eV/ET-) degrees so that 
Pdk ~ 2 P 0 hj Lt. In this work, we do not adopt these approximations. Here, the implicit energy-dependence of the 
UAS rate in Eq. 0 arises from the energy-dependences of Wth, Ai,, and r, as well as from the differential flux Pi,^. 
We present results for the full nested integrals of Eq. 0 - 


B. Horizontal Air-Showers (HAS) 

We now turn to the derivation of the HAS event rate. Neutrino-induced air-showers come in several topologies 0. 
All three neutrino flavors contribute equally to the neutral current (NC) events, but these transfer on average only 
20% of the incident energy to the shower. Furthermore, the NC interaction rate is smaller, about 44%, than the 
charged-current (CC) rate. Among the CC events, the leading muon and tau from incident and i/r, respectively, 
are not visible in the air (unless the tau decays in a “double-bang” event). In the CC process, only 20% of the 
incident energy is transferred to the visible shower. For a zye-initiated CC event, the produced electron contributes 
electromagnetically to the shower, so the full incident energy converts to shower energy. In summary, about one event 
in four (the i^e CC interaction) will transfer 100% of the incident energy to the shower, while three events in four will 
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transfer ~ 20% of the energy. ^ 

To be definite, we assume a r'e CC interaction in what follows. We label the spatial axes as z for vertical upward 
and X and y for the directions tangent to the Earth’s surface; curvature of the Earth’s surface may be neglected 
for HAS. It is useful to consider first a parallel neutrino flux perpendicular to x, incident with a zenith angle 9z, as 
illustrated in the projections of Figs. Inland01 The horizontal air shower probability is 

dH J crJr^Patm(nnt)d^nnt = Tbe 27rACT^^Patm(0) J d COS 0^ J e“^‘“‘/^dzi„t, (13) 

where ri„t is the point of interaction, A = J dx J dy again, and the second expression follows from the first when the 
atmospheric density function 

Patm(2) = /Oatm(O) (14) 

is inserted. We set the scale-height h = 8 km. The factor of f dH = 2 tt J d cos 02 in Eq. (1 1311 rotates the incident 
flux, initially assumed to be parallel and now assumed to be isotropic, over the full sky. 

In principle, we should include the curvature of the Earth’s surface in assessing the vertical height z in p{z) along 
the developing shower. In practice, this is unnecessary as long as the shower length is a small fraction of the Earth’s 
radius, as is the case (we will return to the curvature issue later in the discussion of UAS events). 

The HAS event rate scales linearly with the cross-section cr^^. This is because the absorption probability of the 
neutrino in the atmosphere is negligibly small. The natural scales of atmospheric column density are the vertical 
density 


i?,,(HAS) = J 


dzpatm(z) = /ipatm(O) = 1030 g/cm^ 


(15) 


and the horizontal density 


dhor — 


pOO I nCX) 

/ da;patm(-2; = -v/d?® +a;2 - i?®) « dvert / = \/7ri?©/2 dvert = 36 dv 

do ’ do 


/o ' Jo 

In terms of the latter, the neutrino absorption probability in the atmosphere is 

d 


P(v — air absorption) = Na dhor 


dho 


— 2 X 10 ^ (T 31 


dho 


(16) 


(17) 


where d < dhor is the column density of the neutrino’s trajectory in the atmosphere. Thus, for ^ 10 ^®cm^, 
atmospheric absorption is negligible even for horizontal neutrinos, and so the neutrino interaction rate scales linearly 
with crCC_ 

Further restrictions on the integration variables result from further assumptions for detector efficiencies. Let us 
assume that the air shower must originate in the detector FOV of area A. Then the straightforward integration of 
Eq. gives 

R.^HAS) = 2ttAF,^ da)?^patm(0). (18) 

The value /9atm(0) = 0.62 x IO^'^ctsi sets the scale for the interaction probability in the atmosphere per 

incident neutrino. The resulting value of the acceptance ® is Acc = Rv^{HAS)/Fy^ = 27r A/ipatm(O) = 
3.9(731 ( io 4 ^m^ ) km^sr. This value suggests that wide-angle, large-area detectors exceeding 10"^km^sr, and cos- 
mic neutrino fluxes exceeding 1/km^sryr, are needed for event collection. Put another way, full sky coverage of an 
air mass of ^ 10^ km^ x h p{0) ^ teraton is required. 


^ If the incident neutrino spectrum is falling as a power, then at fixed energy the lAe CC events dominate the total rate. 

® One may also write the acceptance as 27r A (h/Ai/), where Ai/“^ = patm(O) is the neutrino MFP. This expression is the ^ h limit 
of Acc = 277 A (1 — ). In this latter form, one sees the acceptance saturating its geometric value of 27r A in the strong cross-section 

limit. 
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III. CONSTRAINTS FROM DEVELOPMENT AND IDENTIFICATION OF SHOWERS 

In this section we address conditions for the showers to be observable. First of all, shower detection will require that 
within the FOV, the length of the shower track projected on a plane tangent to the Earth’s surface (as would be seen 
from far above or far below) exceeds some minimum length, /min- Space-based observatories are far above the Earth, 
and so view the atmosphere as a two-dimensional plane. For ground-based observatories, tangential projections may 
not be the optimum way to describe the FOV constraint, but we use it as a guide. 

In addition to the projected length constraint, there are three “shower-development” constraints to be applied to 
the events. A minimum column density, dmini beyond the point of shower initiation is required for the shower to 
develop in brightness. On the other hand, after a maximum column density, dmax, the shower particles are below 
threshold for further excitation of the N 2 molecules which provide the observable fluorescence signal. We therefore 
terminate showers at dmax, which implies a finite length for the visible shower. 

While the requirements of minimum projected length, and minimum and maximum shower column-densities are 
correlated, no two of them implies the third. Some reflection on the 9 and z dependences of the varying densities and 
projected lengths reveals that this is so. 

Finally, the fluorescent emission per unit length of the shower will decline exponentially with the air density at 
altitude. At z = 2/i, the fluorescent emission is down to e~^ = 14% of that at sea level. At z = 3/i(4/i), it is 
down to 5% (2%) of that at sea level. Atmospheric absorption of the emitted fluorescence also affects the signal. 
This absorption is thought to scale roughly as the atmospheric density, up to about 20 km [^. Thus, it turns out 
that the fluorescence signal could roughly be taken as constant between zero and 20 km. Accordingly, we will take 
-Zthin = 3/i = 24 km as the “too-thin” altitude beyond which the signal becomes imperceptible. 

Thus, there are four constraints that render the shower observable. These are the ^min, dmin, dmax, and “too-thin” 
(or Zthin) conditions. The values which we choose for these parameters are listed in Table UTl The choice for the Zthin 
value was discussed and motivated above. The dmin and d^ax choices are inferred from the observed longitudinal 
development profiles of ultralngh-energy cosmic ray showers (the famous Fly’s Eye event at energy 3 x 10^° eV 
provides a splendid example Q). Showers at 300-400 g/cm^ of column density (also called “atmospheric depth” 
or “slant depth”) comprise tens of billions of electrons, with a brightness roughly 10% of shower maximum. The 
electrons in showers at ^ 1200 g/cm^ are ranging out, reducing significantly the shower brightness. We assign a 
relatively small value to Imin to maximize the observable event rate. For the EUSO experiment, each pixel is a map of 
a square kilometer of the Earth’s surface |3iT|| . Thus, an Imin of 10 km corresponds to a signal in ten contiguous pixels. 
The background for ten contiguous pixels should be small. With ten pixels, the angular reconstruction of the event 
direction is roughly 1/10 radian (^ 5 degrees). With a cloud layer a smaller value for event triggering may be 
needed. With an ly^in of 5 km, the signal/noise should still be acceptable. For Zmin = 5 km, the angular reconstruction 
is reduced to 1/5 radian (~ 10 degrees), though. Table [171 contains a summary of the values which we choose for the 
four shower-development parameters Zthin, dmin, dmax, and Imin- 

A. Effective area 

Let us describe the projected length Imin constraint, and the dmin and dmax constraints in the general case. Consider 
a detector with a FOV characterized by a radius i?FOV) be., (see Fig. lU 

+ y'^ ^ Rpov ■ ( 19 ) 

A shower produced with initial coordinates (x, y, z) cannot have a visible projected length larger than the chord length 
y + Vd^FOV ~ 11^® FOV. However, the projected length may be smaller, for three reasons. The first is that the 

shower must develop and brighten before becoming visible. This requires traversing the column density dmin. The 
second reason is that the shower may hit the ground (attain the “too-thin” altitude) before reaching the far boundary 
of the FOV in the case of HAS (UAS). And thirdly, the shower may extinguish before reaching the far boundary of 
the FOV. Extinction occurs when the traversed column density attains the value dmax- For HAS (UAS), we label 
the upper altitude where the shower becomes visible (extinguishes or strikes the “too-thin” altitude) as zjj, and the 
lower altitude where the shower extinguishes or strikes the ground (becomes visible) as z^. We will use the HAS or 
UAS label on zu and zz, to distinguish between these altitudes in both cases. The altitudes zb and ze, to be defined 
shortly, will also carry a HAS or UAS label. For brevity, we will sometimes omit the HAS and UAS labels when it is 
clear from the context which label applies. 

From the above shower-development considerations, the visible shower length projected on the Earth’s surface is 
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{zu — zl) tan0. Collecting the remarks above, the projected air-shower length within the FOV of the detector is then 

l-po\ = min |[z (7 - Zl] tan 6*, y -I- \J i?|ov “ | • (20) 

We will discuss the maximum and minimum altitude values zu and zl for each type of shower in the following 
subsections. 

For the projected length in the FOV to exceed some minimum length, Zmin, we infer from Eq. the set of 
conditions 


2/ -l- y Ryow 


( 21 ) 


[zu - Zl] tan0 > • 

After algebraic manipulation of constraints m and the area integral J dx J dy is easily done, yielding 


2 ( arcsinry-77V^l-772 

— TT iL-u 


^FOV 


-K12 


where 




4i?FOV 


( 22 ) 


(23) 


(24) 


is nearly one in a large-area experiment. This is then the constraint-modified meaning of A. Since the arguments 
leading to it apply equally to the HAS and UAS geometry, the result in Eq. applies in both rates, Eqs. 0 and 
m- We note that for Z„iin ^ 2i?FOV) the expression in parenthesis in is nearly unity, with an expansion 


1 - 


/ ^min \ , Q f \ 

\ 2i?FOV / \ 2i?FOV J 


( 25 ) 


Thus, with Zmin ^ 2i?FOVj the constrained area is just the geometric area. 
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B. Four constraints for HAS 

In this subsection we develop the ^minj dmim dmaxi and “too-thin” (or zthin) constraints for HAS. 
The “too-thin” constraint at high altitude is simple. It effectively implies (wint = -^int/ cosO^) 


.S^ground ^ ^int ^ ^thin ? 


( 26 ) 


where Zground = 0 labels the altitude of the ground. We will retain the symbol Zground even though it is zero in the 
present context. Retaining Aground will be useful for substitutions in later sections where we include a cloud layer. 
The projected track length condition, given in Eq. lEU, is 


[zc/(HAS) - zl(HAS)] tan0, > Un • (27) 

Now we consider the calculations of the maximum and minimum visible-shower altitudes, ^[/(HAS) and Zi(HAS), 
respectively, which enter this formula. The column densities of the HAS showers evolved from initial altitude Zint to 
the lower altitude z are given by 

d(HAS; z) = Patm(^) = ^ . (28) 

Jj, cos &z cos &z ^ ' 

The evaluated integral uses the exponential decrease of the atmospheric density with increasing altitude and the 
definition dvert = ZiPatm(O). Equating c?(HAS; z) to dmin, we get the altitude where the shower first becomes visibly 
bright. We call this altitude zb(HAS). Positivity of the integrand ensures that zb < z-mt for any atmospheric density 
profile. Implicitly, zb(HAS) is given by 


g-ZB(HAS)/?l _ g-Zint/^i 


COS Oz ■ 


(29) 


Explicit solutions for zb(HAS) as a function of Zint and vice versa are 


ZB (HAS) = -h\n 



; Zint =-din 

g-ZB(HAS)/?. "^““cosd. 

uvert 


Covert 


(30) 


A shower becomes visible above the ground level only if zb > Aground- From Eq. (I29II . the condition zb > Aground 
implies 


Zint > -h\n 


^vert 


(31) 


Although this condition ensures a visible shower above the ground, it allows the visible length to be arbitrarily small. 
The requirement of a visible projected length in excess of Imin will lead to a stronger constraint, presented below in 
Eq. iHIHll . However, Eq. m is useful in that it implies an absolute limit on the shower direction. Since Zjnt is less 
than zthin by construction (Eq. (l^ b the limit is 

cosd, < ^ ^ . (32) 

^min ^ ^min 

There is no absolute restriction on angle from the dmin constraint if dmin is less than dvert = 1030 g/cm^, since then 
even a vertical shower traverses enough column density to brighten. We will take 400g/cm as our standard value for 

dmin ■ 

The calculation of zb(HAS) proceeds analogously to the calculation of zb(HAS). Setting the column density 
d(HAS; z) equal to dmax, we get the altitude where the shower extinguishes. We call this altitude zb(HAS). Implicitly, 
zb(HAS) is given by 


g-z_E(HAS)//t _ g-zi„t//t 


COS Oz . 


(33) 


The LHS of Eq. (I33II is maximized by setting zb(HAS) equal to Zground and Zint to Zthin- If the RHS exceeds this 
maximum value, then for any Zint the total column density remains less than dmaxi the shower does not extinguish, 
and there is effectively no dmax constraint. So, for more-vertical showers obeying 


cosd^ > /'e-%»und//* _ g-z=thi„//«A = 

<^max ' 


( 34 ) 
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the shower hits the ground before extinction and thus, we have 2 :l(HAS) = ^ground- 

On the other hand, when cos^^ < cos 02, then whether or not the shower extinguishes before striking the ground 
depends on the height in the atmosphere at which the shower originated, i.e., on Zint- Solving Eq. (ESI explicitly for 
0 _e(HAS) as a function of Zint and vice versa, one gets 


ZB (HAS) = -h\n 


+ j”"" cos 02 ] 

; Zint = -/lln 

g-.t.(HAS)//r COS 02 

Overt 


Overt 


= z(HAS). 


(35) 


The shower strikes the ground if zb < -Zground, and extinguishes if zb > z:ground- The critical value is zb = Aground- 
At this critical value, (ISSl) gives 


■Zint = -h\n 


^COS02 

dvert 


= z(HAS). 


(36) 


Thus we have two cases: for Zint < ,z(HAS), the shower strikes the ground and the minimum altitude is zb = Zground; 
while for Zint > z(HAS) (which implies cos02 < cos02 because Zint < -Zthin by construction), the shower extinguishes 
above the ground and zb = zb(HAS). The latter case corresponds to Eq. ll^ having a real-valued solution in the 
physical interval [Zground,-Zthin]) whereas the former case corresponds to no such solution for Eq. II35II . The high 
altitude Z( 7 (HAS) is where the shower begins its visible track length. Accordingly, we set Z( 7 (HAS) = Zb(HAS) in the 
^min constraint for both cases, where zb(HAS) is given in Eq. iHinil . 

For Zint > .z(HAS), the shower extinguishes and we substitute Z( 7 (HAS) = zb(HAS) and zb(HAS) = zb(HAS), 
given in Eqs and (I35II . respectively, into Eq. itTTIl . After a bit of algebra, one finds that the resulting Zniin constraint 
can be expressed as 


-Zint > -h\n 


cos 02 

/ dmax - dmine'*t“»^ \ 

^vert 

1 ^min f 

\ ^ h tan Oz _ j 


(37) 


Real-values of Zint in the interval [z(HAS), Zthin] which satisfy this equation, if any, satisfy all four constraints for 
cos 02 < cos 02 , and so contribute to the integral for the observable event rate. 

For the other case, where Zint < .z(HAS), the shower strikes the ground. We have zb(HAS) = zb(HAS), the latter 
given in Eq. (13011 . and zb(HAS) = Zground- Inputing these expressions into Eq. II27II . one finds an explicit expression 
for the Zniin constraint, 


Zint > -/iln 


g (taHTj: +'^ground)/^ 


drniTl 


COS 02 


(38) 


This constraint ensures a visible projected length exceeding Imin- It replaces the constraint of Eq. EU, which ensured 
only a nonzero visible track length. Of course, in the limit /min = 0, the two constraints are identical. Real-values 
of Zint in the interval [zground, Zthin] for cos02 > cos 02, or in the interval [zground,-z(HAS)] for cos02 < cos 02, which 
satisfy this equation, if any, satisfy all four constraints and so contribute to the integral for the observable event rate. 

To summarize the HAS rate formulas in the absence of cloud cover, we have the general rate equation, Eq. 
with the area given in Eq. (E3, and the “too-thin” constraint in Eq. (P|l. There are two alternate ways to express 
the dmin, dmax) and /min Constraints. The first way is to define the boundaries of (zint, 0z)-integration physically but 
implicitly. This is done with Eq. m implementing the /min constraint, where zj/(HAS) = zb(HAS) implements 
the dmin constraint with zb(HAS) given in Eq. (13011 . and zb(HAS) = max{zground, zb(HAS)} implements the dmax 
constraint with zb(HAS) given in Eq. (I35II . If Eq. II35II has no real-valued solutions in the interval [zground 7 Zthin], then 
Zb(HAS) = Zground- The value of Zground is zero (in the absence of clouds). 

The alternative way to express the dmin, dmax, and /min constraints is to solve the constraints of the first approach 
for explicit boundaries on the (zi„t, 02)-integration. The results of this approach bifurcate, depending on whether the 
shower extinguishes, or the shower strikes the ground. For the case where the shower extinguishes, the boundaries are 
given by Zint > -z(HAS), with z(HAS) defined in Eq. (1^ . and by Eq. E3- For the case where the shower strikes the 
ground, the boundaries are given by Zint < -z(HAS), and by Eq. (I38II . A priori, there is no guarantee that Eqs. E3 
and have real-valued solutions in the physical region of Zint- 


C. Four Constraints for UAS 

The calculation of the /min, dmin, dmax, and “too-thin” (or Zthin) constraints for UAS events proceeds analogously 
to the calculation for HAS events. The “too-thin” altitude constraint is again 


Zground ^ Zdk ^ Zthin ■ 


( 39 ) 
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At sea level, Zground = 0. The projected track length condition for UAS events, analogous to Eq. (ETIl for HAS events, 
is 


[z; 7 (UAS) - 2 ;l(UAS)] tan6»„ > /min ■ (40) 

The values of zc/(UAS) and 2 ;l(UAS) differ from Z[/(HAS) and zl(HAS). To calculate them, we turn to calculations 
of UAS column densities, given by 


d{VAS;z)= j 

J ZA 


dz 

cos 6, 


cos On V ^ 


Setting this equal to dmin defines implicitly the brightness altitude zb(UAS): 


_ g-ZB (UAS)//^ = ^ cos, 


(41) 


(42) 


Solving this equation explicitly for zb(UAS) as a function of Zdk and vice versa, one gets 


2 :b(UAS) = —din 



; Zdk = -Iiln 

[j““cos0„ + e-^®(UAS)/J 

avert 


avert 


(43) 


We require that zb(UAS) < Zthin; otherwise, the shower invisibly disappears into thin air. From Eq. (I43II . the 
condition 2 :b(UAS) < Zthin can be written 


Zdk < —/iln 


dmiTi 


■ COS 


On + e 


(44) 


This condition ensures visibility of the shower, but with a visible length arbitrarily small. The requirement of a visible 
projected length in excess of /min will lead to a stronger constraint, presented in Eq. below. Positivity of Zdk and 
Eq. 1331) lead to the same angular constraint for cosdz as was found for HAS’s cosd„ in Eq. 1221. 

Setting d(UAS; z) equal to dmax, we get the high altitude 2 :_e(UAS) where the UAS shower extinguishes. Implicitly, 
this highest visible altitude is given by 


g-zdk//« _ g-ZE(UAS)//« ^ ^ 

dvert 


(45) 


The LHS of Eq. is maximized by setting ^dk equal to Aground and z_e(UAS) to 2;thin- If the RHS exceeds this 
maximum value, then for any Zdk the total column density remains less than c/max, the shower does not extinguish, 
and there is effectively no dmax constraint. So, for more-vertical showers obeying 


(46) 


we have zt/(UAS) = zthin- 

On the other hand, when cosd„ < cos0„, then whether or not the shower extinguishes before reaching the “too¬ 
thin’’ boundary depends on the height in the atmosphere at which the shower originated, i.e., on 2;dk- Solving Eq. 1321 ) 
explicitly for 2 £;(UAS) as a function of Zdk and vice versa, one gets 


Zi5(UAS) = -hln 



; Zdk = -^ In 

'^““cos0„ + e-"®(UAS)A 

(2vert 


_ Covert 


(47) 


The shower extinguishes if z_e(UAS) < zthin, and hits the “too-thin” boundary if Z£;(UAS) > Zthin- The critical value 
is 2 :£;(UAS) = ^thin- Inputting this critical ze into Eq. lUTI) . one finds a critical value for the decay altitude 


Zdk = -hln 


■ cos 


0n + 


= z(UAS). 


(48) 


For Zdk < z(UAS), the shower attains the “length” c?max and extinguishes, whereas for Zdk > z(UAS), the shower 
reaches the “too-thin” boundary zthin without extinction. 
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For Zdk < -?(UAS), the shower extinguishes and so zj/(UAS) = 2 :£;(UAS). Substituting this (Eq. (ETIl i and 
zl(UAS) = zb(UAS) from Eq. (Iddll into the ^min constraint Eq. (EUI) leads to an explicit expression for the ^min 
constraint: 


Zdk > —hln 


cos 9n 

(d ph tS; - d ■ \ 

' ^max ^ ^min \ 

^vert 

1 imin I 

\ g tan Bjx _ \ ) 


= z< 


(49) 


This limit, like the analogous one for HAS in Eq. m, forces the shower initiation to occur at a higher altitude where 
the air is thinner, and therefore, for fixed dmaxi the shower and its projection are longer. There is no dmin constraint 
for a shower that saturates dmax- 

For Zdk > z(UAS), the shower reaches zthin without extinction. Therefore, there is no dmax constraint and we set 
2 ;c/(UAS) = 2 :thin- Substituting this and 2 :l(UAS) = zb(UAS) into the Zmin constraint Eq. (BHll . we find the following 
explicit expression for the Zmin constraint: 


Zdk < —din 


+ 


■ cos 


= 


(50) 


For this class of showers which reach the “too-thin” boundary, this upper bound on Zdk ensures a visible projected 
length exceeding Zmin. It supersedes the constraint of Eq. (EH), which ensured only a visible track of non-zero length. 
Of course, in the limit /min = 0, the two conditions are identical. 

The essence of the /min constraint is that the shower must have sufficient normal angle to attain a minimum 
horizontal projection. Thus, there is a critical angle 0“'* for which the conditions of Eqs. II49I) and collapse to 
z< = z(UAS) = z>. For normal angles smaller than 0“'*, i.e., for cos 6 („ > cosd™*, there are no observable events. 
Setting z< = i(UAS) = z>, one finds that this critical angle is given implicitly by 


zthin = -hln 


dmax dmin 


cosd™* 

^min 

g h tan ^ 


(51) 


This critical angle encapsulates a relatively weak constraint in the cloudless case. However, it will become a strong 
constraint when we consider cloudy skies. 

To summarize the HAS rate formulas in the absence of cloud cover, we have the general rate equation, Eq. 0 , 
with inputs from Eqs. (0), (©, ll51)-(tni) and (EH, and the “too-thin” constraint in Eq. EH- As was the case for 
the HAS events, there are two alternative ways to express the dmin, dmax, and /min constraints for UAS events. 
The physical but implicit approach bounds the (zdk, dn)-iiitegration with Eq. 11401) implementing the /min constraint, 
where zj/(UAS) = min{zthin)-Z£;(UAS)} implements the dmax constraint with Z£:(UAS) given by Eq. (1471) . and where 
zl(UAS) = Zb (UAS) implements the dmin constraint with zb(UAS) given in Eq. (ESI). 

Alternatively, one can solve the constraints of the first approach explicitly for the (zdk, dn)“iiitegration boundaries. 
As with the HAS events, the results again bifurcate, depending on whether zb(UAS) < Zthin (the shower extinguishes), 
or zb(UAS) > Zthin (the shower runs out of air). For the case where the shower extinguishes, the boundaries are given 
by Zdk < z(UAS) with the latter quantity defined in Eq. (HHl) . and by Eq. (EH- For the case where the shower strikes 
the “too-thin” altitude, the boundaries are given by Zdk > z(HAS) and by Eq. (03) • 


D. Four Constraints for UAS, Including Earth’s Curvature 

So far we have treated z as the vertical height above a flat Earth. As remarked in the section on HAS rates, this 
is a valid approximation as long as the trajectory in the atmosphere is small relative to the Earth’s radius Rq. Such 
is the case with HAS events. However, at 10^° eV, the decay MFP for a tau is nearly 5000 km, comparable to i?® 
(6371 km). At extreme energies, curvature effects cannot be neglected for UAS events. Specifically, the error made 
in the vertical height, as a function of the atmospheric path length w and angle dhoD is (neglecting terms of order 

OiwVRl)): 


w'^ COS^ dhor 

2 i ?0 1 -|- sin dhor 2 i ?0 


(52) 


The far RHS expression, ui^/ 2 i ?0 = 78(r(;/10^km)^ km, is saturated for (near) horizontal events. The error 5z in 
the height of the decaying tau can be considerable: ~ 2000 km for nearly horizontal events with Er ~ 10^° eV, 
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FIG. 5: Illustration of the dependence of shower variables on the Earth’s curvature. For clarity, the various parts are not 
drawn to proportion. 


and ~ 20 km for nearly horizontal events with Et ~ 10^® eV. So our neglect of Earth’s curvature means that we 
underestimate the height of the shower, and so overestimate the air-density for shower development. The height 
underestimate will erroneously reduce (increase) the event rate as viewed from space (ground) when we introduce 
clouds. The density overestimate will erroneously enhance shower development. 

The curved geometry is shown in Fig. The net effects of curvature are twofold. First, for a given pathlength 
w and trajectory angle 9^ at emergence from the Earth, the curvature-corrected altitude, which we label as z\ is 
increased. Second, the angle of the shower with respect to a plane tangent to the Earth directly below, which we 
label as 9'^^^ and 9'^ for the horizontal and nadir angles, respectively, are rotated relative to the comparable emergence 
angles (again, for a given w and 9n)- Note that the primed variables are the altitude (z') and angle {9') seen by a 
detector. The unprimed variables describe the shower’s prehistory. 

We now list the geometric relations which we need. From applying the Pythagorean theorem to the right triangle 
in Fig. 13 we get z' in terms of w and 0hor: 




(53) 


Applying the Law of Sines to the same right triangle, we get 

sin 9n cos 0hor 


sin 9'„ = cos = 


l + ir l + ir 


with z'{w,9hoT) given in Eq. lO. From this comes 


cos 9'„ = sin = 


sin ^hor 




1 + — 


/ 5 


and 


cot 6'^ = tan = tan 6*hor + 


w 


cos ^hor 


(54) 


(55) 


(56) 


This last expression shows clearly that for nearly tangent (“Earth-skimming”) events, the observed angle (^^qj.) 
increased from the emergent angle by the term « w/Rq. 
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In particular, the altitude and the angle at the decay point of the tau are obtained by setting w equal to u>dk in 
Eqs. (ICT - lIKHll . For example, the decay altitude z'^y. is given by 



(57) 


After the tau decays, the resulting shower has a length which is short on the scale of i?®. Consequently, we may 
ignore the Earth’s curvature from Wdk onward, and use and 0'^ as the shower angles. 

The development of the shower constraints with curvature parallels that for UAS events without curvature, but 
with rather than {z,0n) parameterizing the shower altitude and angle. 

The “too-thin” shower-constraint now becomes 


'^dk(^dk; ^hor) ^ 

Using Eq. 113, this constraint can be cast as a constraint on the integration variables Wdk and 0hor (equivalently, L). 
The result is 


^dk ^ ^dk (^hor; -^thin) ? 


(59) 


with 


^dk 

-^0 


I sin ^hor 


/ -2^thin \ 

\WJ 


2 + 


-^thin 


-Re / 




sin ^hor • 


(60) 


This constraint requires the decay to occur at an altitude below zthin, but does not require a minimum of shower 
development or a minimum of projected length. The shower development and /min constraints to come will be stronger. 
However, one useful feature of this constraint is that it defines the maximum decay distance available to a tau in our 
curved atmosphere. Taking 0hor = 0 to maximize the available decay length, one finds 


iw“k< w)(j,“(6>hor = 0) = 0Zthin (2i?e + -^thin) ^ 550 ykm. (61) 

The mean decay length for a tau is already 490 km at 10^® eV, and grows linearly with its energy. Thus, it is clear 
from Eq. EJ that the Earth’s curvature significantly reduces the UAS rate at energies at and above 10^® eV. It is 
also clear from (EH) that since <C R®, leading order expansions in the ratio w/Rq are valid. For example, to 

very good approximations, Eqs. lEKll and E3 can be written as 

sin = sin 6'hor + ^ , (62) 

-^0 

and 

2 

4k ~ '^dk sin 6»hor + 7^ ■ (63) 

Zn^ 

The first term on the RHS is just Zdk for a “flat Earth”. So we learn that the replacement of Zdk with z'^^ in Eq. 16.‘lll . 
like that of 0hoi (at small angle) with in Eq. ESI), is a simple translation. 

The constraints from shower development require explicit expressions for the shower column density d{wdk, 0n) and 
shower length. The column density in the observer’s primed coordinates is 


d{lJAS;z')= f dwp{z')=f -^^p(z') 

Jz' cos 0'„ 


cos0(j V 



(64) 


with the variables 0'^ and z'jj. depending on just Wdk and 0n- Changing z ^ z' and 0 ^ 0' in Eqs. ij^ and lUKll defines 
2 :^(UAS) and ^^(UAS), respectively. The explicit expressions for 2 ;^(UAS) and ^^(UAS) are given by applying z ^ z' 
and 0n —> 0'n to Eqs. 63 and lITTI) . respectively. The results are 


4 (UAS) = -h\n 



dmiTi 


COS 0' 


(65) 
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and 


z^(UAS) =-/iln 



cos6»; 


( 66 ) 


The Irnin Constraint is 


\Zu tan > ^min 7 


(67) 


with 


z[j = min{zB(UAS), Zthin} , and = Zs(UAS). 


( 68 ) 


An algorithm has emerged for including Earth-curvature in our prior calculations. We simply replace the unprimed 
variables Zdk, On, zb(UAS), z_e(UAS) with 9', Z 3 (UAS), z^(UAS), where z'^^^{zdy,,9n) and 6'^{zdk,9n) are given 
in Eqs. and (UAS) and z^(UAS) are given in Eqs. (I65II and ®. Note that the parameters Zthin, 2 :ground, 

and soon-to-be-introduced Zdoud are never primed, for they define layers concentric with the spherical Earth. 

At this point we can summarize very easily the UAS constraints including the Earth’s curvature (in the absence 
of cloud cover). We have again the general rate equation, Eq. ®, with inputs from Eqs. 0 , 0 , ©-lO and 
(HI. The “too-thin” constraint is given by Eq. (EHli, or equivalently, by Eqs. and (EOll. The Imin constraint is 

given by Eqs. lIHTll and ra . with (UAS) in Eq. lUIKll implementing the dmin constraint and z^(UAS) in Eq. 16611 
implementing the dmax constraint. Unprimed variables are obtained from primed variables via Eqs. , or to 

a very good approximation via the simplified Eqs. In any case, we have used the latter for our numerical 

computations. 

As in the “flat Earth” calculation, more algebra can be done when Eq. EHl is substituted into Eq. The result is 
Eq. lUnil for the l^in constraint if z^(UAS) < z'(UAS), and Eq. (IKHIl if z^(UAS) > z'(UAS), where (zdk, 9n) (z:'ij., 0)^) 

in Eqs. 103 and and z'(UAS) is defined by priming Eq. I48II : 


z'(UAS) = -h\n 


■ cos 9' +e "’“■in/'t 


(69) 


However, when these explicit constraints are expressed in terms of the unprimed variables, no easy separation of Zdk 
and 9n appears to be possible. Consequently, it seems best to treat the constraints as non-linear relations in the 
(zdk, 0)-integration space. This is exactly what we do. 

As a check on our work, we have regained the “flat-Earth” results for UAS from the curved-Earth formalism, by 
taking the Earth’s radius to be very large in the constraint equations (but keeping physical in Eq. 0 and the 
L = 2 i ?0 cos 9n relation). 


IV. SHOWER RATES WITH CLOUDS 

The presence of clouds, their distribution, altitude and optical depth would obviously affect the observed event 
rates. Hence, when we consider the effect of a cloud layer on the observable event rate, the resulting constraints 
become more complicated. In addition, the constraints come to depend on whether the detector is above or below the 
clouds, i.e., on whether the detector is space-based or ground-based. We will model the cloud layer in a very simplified 
way as an infinitely thin layer, but with infinite optical depth and we will assign Zdoud to be the height of the relevant 
cloud boundary. However, we will not take into account how the cloud presence could affect the reconstructed shower 
geometry and energy [^. For a space-based detector, the only visible air showers are those at z > Zdoud, while for 
a ground-based detector only showers at z < Zdoud can be seen. Let us name the four possible event-detection types, 
UAS and HAS as seen from space (S) or from the ground (G), as UASGiS, UAS( 8 )G, HAS( 8 )S, and HAS( 8 )G, in obvious 
notation. Our calculation is partitioned into four parts, corresponding to these four event types. We will see that 
UAS®G and HASoS event types are easily calculated with only simple modifications of our prior, cloudless formulas. 
However, the UASoS and HASoG event types require more care, since the observed shower may have its origination 
above or below the cloud boundary for these cases. 

In our simplified model, the actual development of the shower does not depend on the presence or absence of clouds. 
Thus, the expressions for zb and ze, determined by the dmin and dmax constraints, are unchanged. However, the 
visible projected length of the shower certainly depends on the presence or absence of clouds. 

In this section we do not include the effect of the Earth’s curvature. We have seen that Earth curvature does 
not affect the HAS calculation. In the next section, m we include the Earth’s curvature along with clouds in the 
calculations of the UAS rates. 
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A. UAS in Ground-Based Detectors 

Looking upward with a UAS(8)G detector, the observable atmosphere is bounded from above by the cloud layer. 
In the absence of clouds, the observable atmosphere was bounded from above by Zthin- Thus, the prior, cloudless 
calculation applies to the cloudy atmosphere if we just reset the “too-thin” height Zthin to the cloud boundary ^cioud- 
One feature of this replacement is that the angular constraint on UAS events, traceable back to Eq. (EH, becomes 

cos6i„ < (l - = cose*Q (70) 

^min ^ ' 

If Zcioud < —/iln[((ivert — dmin)/dvert] > then COS 0 q is less than 1, presenting a real constraint on the shower direction. 
With dmin = 400g/cm^, cos^q < 1 occurs for Zdoud < 3.9 km. Thus, for cloud boundaries below this value, near¬ 
vertical showers do not satisfy the dmin condition. For example, with a cumulus cloud layer at Zcioud = 2 km, only 
shower angles > 55° are allowed by the d„iin constraint. For common extrapolations of the neutrino-nucleon 
cross-section to very high energies, only “Earth-skimming” neutrinos are expected to emerge from the Earth at very 
high energies. For these “Earth-skimming” neutrinos, clouds must be very low to affect the rate. 

We rewrite the relevant UAS constraint equations with zthin ^ -^cioud to include the cloud boundary. The “too-thin” 
constraint in Eq. EH becomes 


0 ^ -^dk ^ -^^cloud ■ 


The ^min constraint Eq. EH is replaced with 


(71) 


[min{z£;(UAS),Zcioud} - 2 :b(UAS)] tand„ ^ ^min 5 (72) 

with Z£;(UAS) and zb(UAS) as before, given in Eqs. (I47II and EH, respectively. This concludes the implicit calculation 
of the integration boundaries. 

Explicit boundaries are obtained by substituting Eqs. m and EH into EH- The result is that the Zmin constraint 
is given by Eq. ll^ when Zdk < z(UAS 0 G), and by (derivative from Eq. (1^ via Zthin ^ -Zcioud) 


Zdk < —hln 


■ cos 9r, + ^cioud)/?!. 


(73) 


when Zdk > z(UAS 0 G); the critical altitude z(UAS 0 G) (derivative from Eq. (1^ via zthin ^ -Zcioud) is given by 


z(UAS® G) = -din 


g ^cloud/^ I 


COS 9n 


(74) 


As with the no-clouds case, these /min constraints require a sufficiently large normal angle so that the horizontal 
projection of the shower is visible. The critical angle is obtained from Eq. EH with the substitution zthin by ^doud- 
The resulting equation is 


■^cioud — h In 


^max 


cos6»^"^^ 

^min 

g h tan ^ 


(75) 


The meaning is that given a cloud layer at Zcioud, there are no visible events for cos0„ > cos0^^*. Since Zcioud is a 
monotonically increasing function of cos0“'*, this result may be stated in a different way: there are no visible events at 
cos 0n > cos 0“'* if there is a cloud layer lower than that of Eq. EH- There are no visible events at all if cos0™* = 0. 
From Eq. EH, this occurs for the critical cloud altitude 

dmax djnin ^ 

—d- r~ ■ 

^vert ^min_ 

Thus, clouds completely obscure the detector if Zcioud < -Zcfoud- 

To summarize the formulas giving the UAS(8)G events even in the presence of cloud cover, the general rate equation 
is Eq. El, with inputs from Eqs. (EJ, (0, ®-EH 12H, the “too-thin” constraint is given by Eq. EH, the /min 
constraint by Eq. EH, with Z£:(UAS) and Zb(UAS) unchanged from their cloudless expressions, Eqs. II47II and EH, 
respectively. Alternatively, explicit /min constraints are available in Eq. (I49II for Zdk < .2 (UAS 0 G), and in (I7dll for 
-sidk > z(UAS 0 G), with z(UAS 0 G) defined in Eq. (17^ . These twin constraints reflect the two possible outcomes of 
min{zs(UAS), Zcioud} in Eq. EH- 


crit _ _ t 1 
•^cloud ^ 
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B. HAS in Space-Based Detectors 

Looking downward with a HAS(8)S detector, the visible atmosphere is bounded below by the cloud layer. In the 
absence of clouds, it is bounded below by ^ground- Thus, the cloudless calculation applies when Aground is reset to 
^cioud- With this type of substitution in mind, we retained the symbol Aground in our prior cloudless formulas, even 
though its value was zero there. For example, for the LIASES events, the angular constraint of Eq. (EH becomes 

cos ^ ^ g* (77) 

^min ■' 

This constraint ensures that the shower brightens sufficiently above the clouds to become visible. If Zdoud > 
—/iln(di„in/dvert+ then cosdg < 1, being a real constraint on the shower direction. With dmin = 400g/cm^, 

cos 0g < 1 for Zcioud ^ 6-6 km. Cumulus cloud layers rarely rise to this height, and so there is no angular constraint on 
HAS resulting from cumulus clouds. However, cirrus clouds populate the high atmosphere, and therefore do constrain 
the HAS angle 9z- 

We rewrite the other relevant HAS constraint equations with Zground Aloud to include the cloud boundary. The 
“too-thin” constraint becomes 


Aloud ^ ^int ^ Ahin ■ 


(78) 


The ^niin constraint Eq. EH becomes 

[(zb(HAS) - max{z_E(HAS), Zcioud})] tan^^ ^ ^min 5 (79) 

with zb(HAS) and Z£;(HAS) as before, given in Eqs. UdOB and 118511 . This concludes the implicit calculation of the 
integration boundaries. 

Explicit boundaries are obtained by substituting Eqs. EH and (I85II into EH- The result is that the Zmin constraint 
is given by Eq. (I87II . for Zjd > i(HAS 0 S), and by the following (derived from Eq. (1881) via Zground ^ Aloud) for 
Ant < z(HAS 0 S): 


Zint > -/lln 


- Poland)/h _ “min 


COS 9, 


the critical altitude z(HAS 0 S), derivative from Eq. Il8til) via Zground ^ Aloud, is 


z(HAS(g)S) = -h\n 


g-.ctoudA_^COS0, 

^^vert 


(80) 


(81) 


To summarize the formulas giving the HAS(8>S events even in the presence of a cloud layer, the general rate equation 
is Eq. EH, with the area given in Eq. EH, the “too-thin” constraint is given by Eq. EH, and the In^in constraint by 
Eq. (EH, with Zb (HAS) and zb(HAS) unchanged from their cloudless expressions. Alternatively, explicit expressions 
for the /min constraint are available in Eq. lIHIIll for Zint < z(HAS (g) S), and in Eq. (TTtIi for zint > z(HAS 0 S), with 
z(HAS 0 S) defined in Eq. (1811) . These twin constraints reflect the two possible outcomes of max{zB(HAS), Zdoud} in 

Eq. EH- 


C. UAS in Space-Based Detectors 

Looking downward with a UAS^S detector, the visible atmosphere is bounded below by the cloud layer, but the 
visible UAS may have begun its development above or below the clouds. Thus, the “too-thin” constraint remains 
Eq. (1891) as in the cloudless calculation. The new /min constraint is 

[min{zB(UAS), Zthin} - max{zB(UAS), Zdoud}] tan6»„ > /min (82) 

with zb(UAS) and zb(UAS) as before, given in Eqs. (H7I) and EH, respectively. This concludes the implicit calculation 
of the integration boundaries. The left-hand side of Eq. (EH makes it clear that generation of a visible UAS requires 
Zb(UAS) > Zdoud and zb(UAS) < zthin, and that the clouds are irrelevant when zb(UAS) > Zdoud- 
Explicit boundaries, if desired, are obtained by substituting Eqs. (IT7I) and EH into (1^ . There are 4! = 24 a 
priori orderings of the four parameters in the /min constraint. However, the orderings zb(UAS) > zb(UAS) and 
Ahin > Zdoud are fixed. This leaves 4!/(2 • 2) =6 possible orderings of the parameters. Of these six orderings, one 
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has zb(UAS) > Zthin and another has 2 ;£;(UAS) < Zdoud- These orderings do not produce an observable shower, the 
former showering too late and the latter showering too early. We are left with four relevant orderings: 

(a) Zthin ^ ^ ^ -^cloud; 

(b) Z]^ > '2:thin ^ ^ ^cloud; 

(c) .^thin ^ ^E ^ ^cloud ^ ^B: 

(d) > .2:thiii ^ .2:cloud ^ Zb • 

The first two, (a) and (b), are characterized by zb > -Zcioud, which restricts Zdk according to 


2dk > —ft. In 


g-.c.oud/?> + COSffr, 


(83) 


Ordering (a) characterizes the shower that extinguishes (i.e. Zdk < i(UAS)), while (b) characterizes the shower that 
reaches the “too-thin” air boundary (i.e. Zdk > .z(UAS)). For the two orderings (a) and (b), the clouds do not obscure 
any part of the visible shower, and the /min formulas of the cloudless section Hill 01 apply. The next two orderings, (c) 
and (d), are characterized by Zdoud > zb. For these two orderings, the clouds do obscure part of the visible shower. 

However, it may be that for low values of Zdoud, there are no events satisfying the topologies specified in (c) and (d) 
regardless of whether the cloud layer is actually present at Zdoud- For example, with our canonical value Zthin = 3ft, 
the conditions in (d) require that the shower remain visible over a vertical length at least as long as (3ft — Zdoud)- 
For a small value of Zdoud 7 such a shower will never happen, as showers cannot both begin below Zdoud and survive 
beyond 3ft. Thus, in order to see if there is some critical altitude below which clouds would not suppress the rates 
in UAS(8)S detectors, we seek the conditions for which categories (c) and (d) do not contribute events even in the 
absence of clouds. Under such conditions, the acceptance is just that of the cloudless case. 

For ordering (d), i.e. Z£:(UAS) > Zthin and Zdoud > zb(UAS), the /min constraint presents a restriction on alone. 


[.^thiii Z^doud] tan ^ /min 7 OI" ^doud ^ "^thin Imin/tanftn ■ 


(84) 


For ordering (c), i.e., Zthin > zb(UAS) and Zdoud > zb(UAS), the /min constraint is 

[zb(UAS) - Zdoud] tanft^ > ^min ■ 

Inputting zb(UAS) from Eq. (14711 leads to a restatement of this latter constraint as 


Zdk > —ft In 


■cos6»„ + e (tlne"„+^cioud)A 


(85) 


( 86 ) 


Hence, the altitude on the RHS expresses the minimum altitude above which tau decays contribute events to cate¬ 
gory (c). Since we have argued that the minimum altitude for shower development in category (d) is higher than in (c), 
the RHS also expresses the minimum altitude above which tau decays contribute events to categories (c) and (d). 
Thus, the RHS is the minimum altitude above which tau decays produce showers partially obscured by clouds. 

This constraint on Zdk is analogous to the one in Eq. which sets the minimum value for Zdk in the absence of 
clouds. Since clouds must be more restrictive than no clouds, the RHS of Eq. (IS^ must be larger than z< (RHS of 
Eq. Il4bll l if categories (c) and (d) are to have events. This happens if 


Zdoud ^ ft in 


dmax dniin COS 


- I 


(87) 


This is a necessary condition for events to fall into categories (c) and (d). Consequently, the inequality of opposite 
sign is the sufficient condition for clouds to not obscure the showers. The RHS of Eq. lISTJi is a critical cloud altitude. 

For the typical parameter choices which we consider, the RHS of Eq. ilSTI) . when positive, has a very weak dependence 
on cosOn- In particular, it does not differ much from its value evaluated at cos 9n = 0, which is 

ftmax ^min ft /oo^ 

-3- ] - ■ 

^vert ^min_ 

Thus, we may use Eq. (ISSli as a very good approximation to the RHS of CT . The approximation becomes even better 
as the cross-section becomes larger, for then events come from more horizontal neutrino trajectories. Conveniently, 
the definition of in Eq. is identical to that in Eq. iTTHll . Thus, for all practical purposes we can use the 

same zOf^^j as the critical altitude for space-based and ground-based detectors. ® Thus we have the complementary 


z“‘ud = -h\n 


® We have checked numerically the accuracy of this statement. 
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situation that clouds below ^cioud completely obscure UAS(8)G, but do not affect UAS0S at all. Of course, clouds 
above -2“oud partially obscure both UASoG and UASoS. 

This concludes the explicit construction of the UAS®S constraints for all four allowed orderings of the parameters. 
The final event rate is the sum of the contributions from the four allowed orderings. As a check, we note that in 
the limit Zdoud ^ 0, orderings (c) and (d) no longer contribute, since zb > 0. Thus, we are left with just the two 
orderings (a) and (b) of the cloudless limit. 

To summarize the formulas giving the UAS(8)S events even in the presence of a cloud layer, the general rate 
equation is Eq. 0 , with inputs from Eqs. ©, ®, (ISl-flU). and (EHll . The “too-thin” constraint is given by (TTHll . the 
^min constraint by Eq. (I82II . with 2 :£;(UAS) and zb(UAS) unchanged from their cloudless expressions (I47II and m- 
Explicit solutions for the Zmin constraint, if desired, are given above for the four allowed orderings of the parameters 
Zthin)-s^cioud}- The total event rate is the sum of the four contributions. 


D. HAS in Ground-Based Detectors 

Looking upward with a HAS( 8 )G detector, the visible atmosphere is bounded above by the cloud layer, but the 
visible HAS may have begun its development above or below the clouds. Thus, the “too-thin” constraint remains 
Eq. ll^ as in the cloudless calculation. The Z„iin constraint is 

[min{zB(HAS), Zdoud} - max{z£:(HAS), 0}] ta,n9^ > (89) 

with z_b(HAS) and z_e(HAS) as before, given in Eqs. (TTHll and 13511 . respectively. This concludes the implicit calculation 
of the integration boundaries. Note that the left-hand side of Eq. makes it clear that generation of a visible shower 
requires z_e(HAS) < Zdoud, and that the clouds are irrelevant when zb(HAS) < Zdoud- 

As was the case with UAS 0 S events, there are four orderings that contribute to the HAS( 8 >G events: 

(a) z:c[oud ^ ^ ^ O 5 

(b) Zdoud > ZB > 0 > ze; 

(c) Zb ^ -^doud ^ 0 > Ze , 

(d) ZB > z:doud > ZB > 0. 

The first two, (a) and (b), are characterized by Zdoud > zb- For these two, the clouds do not obscure any part of 
the visible shower, and the Imin formulas of the cloudless section iilll Bl apply (ordering (a) characterizes the shower 
that extinguishes, while (b) characterizes the shower that reaches the ground). The next two orderings, (c) and (d), 
are characterized by zb > ^doud- Here, the clouds do obscure part of the visible shower. For ordering (c), i.e., 
zb(HAS) > Zdoud and 0 > zb(HAS), the Imin constraint presents a restriction on alone. It is 

■ 2 :doud tan Oz ^ ^min • (90) 

For ordering (d), i.e. Zb(HAS) > Zdoud and Zb(HAS) > 0, 

[^doud - z:b(HAS)] tanOz > ^min ■ (91) 


A short calculation leads to a restatement of this latter constraint as 


^int < -hln 


g ( taUal, Zdoud) 


COS 9z 


(92) 


This concludes the explicit construction of the constraints for all four allowed orderings of the parameters. The final 
event rate is the sum of the contributions from the four allowed orderings. As a check, we note that in the limit 
Zdoud —*■ Zthin, orderings (c) and (d) no longer contribute, since zb < Zthin- Thus, we are left with just the two 
orderings (a) and (b) of the cloudless limit. 

To summarize the formulas giving the HAS(8>G events even with a cloud layer present, the general rate equation is 
Eq. n 311 . with the area given in Eq. II23II . The “too-thin” constraint is given by Eq. 12611 . the Imin constraint by Eq. I89II . 
with Zb (HAS) and zb(HAS) unchanged from their cloudless expressions Eqs. I35II and 13011 . Explicit solutions for the 
^min constraint, if desired, are given above for the four allowed orderings of the parameters zb(HAS), Zdoud, zb(HAS), 
and Zground = 0. The total event rate is the sum of the four contributions. 
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FIG. 6: Critical altitude below which clouds would obscure the detector for UAS(8)G and have virtually no effect for UASCiS. 
Above the critical altitude, clouds would partially obscure both UAS(8)G and UAS®S. The label on the abscissa can be thought 
of as an experimental sensitivity to showers. We have taken dvert/h = 1030 gcm~^/8 km. 


E. Remark on UAS with clouds 

The RHS of Eq. ll7?)ll gives the critical altitude below which clouds would completely obscure ground-based detection. 
The RHS of Eq. ll??Hll gives the critical altitude below which clouds would not affect space-based detection. These 
two equations are inverse to each other in meaning, but numerically the critical altitudes on the RHS’s in these two 
equations are identical: 




dmax dmin 


(93) 


Here we have factored the argument of the logarithm into a ratio (dmax — dmin)/^min which is determined by the 
experimental triggers (i.e., by humans and their optics), and the ratio h/d^ert = l/p(0), which is Nature’s gift of our 
atmospheric density at sea level (cf. Eq. (Ilbll l. Numerically, the latter term is 8km/1030gcm“^ = (129gcm“^/km)“^. 
When the argument of the logarithm is < 1, this equation has a positive solution, and so sufficiently low-lying clouds 
will completely obscure UAS0G, but not affect at all UAS(8)S. There will always be a positive solution Zdoud whenever 
(dmax dmin)/^min Is leSS than dvert~ /^(O) — 129 g Cm /km. 

The range of a visible shower at or near 10^° eV is comparable to dvert, and so (dmax — dmin)/dvert is of order 
unity. Typically, the visible length required for shower identification is of order of h = 8 km, so h/lmin is also of order 
unity. Thus, the argument of the logarithm is of order unity. Consequently, whether there can be signihcant cloud 
obscuration in UAS(8>G or no effect in UAS^S will depend critically on an experiment’s choice of shower parameters, 
dmin, dmax, and As relevant examples, with the choices dmin = 400 (300) g/cm^ and dmax = 1200 (1500) g/cm^, 
for Zmin = 10 km the argument of the log is < 1, is 3.8 (0.56) km, and so a cloud layer at a lower altitude 

will completely suppress observation in UAS(8>G detectors and have no effect in UAS(8>S detectors. In contrast, for 
^min = 5 km the argument of the logarithm exceeds unity, there is no Zcioud, and so there is partial rate suppression 
due to clouds at any altitude for both UAS(8)G and UAS(8)S. 

In Fig. El we plot z^d/u^ versus the trigger-parameter combination (dmax — dmin)/^min, over the range 30- 
130 km/gcm“^. As foretold, for values above dvert/h = p(0) = 129km/gcm“^, the solution to Eq. (li^ is negative and 
partial cloud suppression occurs in both UAS(8)G and UAS(8)S. But for (dmax — dmin)/^min less than 129km/gcm“^, 


To be accurate, it is the minimum of the RHS of J87| which gives the critical altitude, as we explained in Section llV (II 
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there is a positive critical altitude delineating total cloud suppression as seen from ground and no cloud suppres¬ 
sion as seen from space, from partial suppression of both. We infer from the figure that an experimental trigger 
(dmax — dmin)/^min exceeding (50,80,100) km/gcm“^ is required for (i) UAS(8)G to avoid complete rate suppression 
from clouds below (7.6,3.8,2.0) km, while suffering partial suppression from clouds above (7.6, 3.8, 2.0) km; (ii) UAS0S 
to have partial rate suppression from clouds above (7.6, 3.8, 2.0) km, while suffering no suppression from clouds below 
(7.6, 3.8, 2.0) km. 

The ground-based result is intuitive, in that as the trigger sensitivity (dmax —dmin)/^min is increased, the experiment 
may tolerate clouds ever closer to the ground. The space-based result is less intuitive. In this case, as the trigger 
sensitivity is increased, the critical cloud altitude above which UAS are partially obscured is again lowered (after 
all, the same is common to UASgG and UAS(8>S). The reason is that with better triggering, a space-based 

experiment is sensitive to more effective atmospheric volume, and so to the presence of lower cloud layers. Put another 
way, with better sensitivity a space-based experiment may see deeper into the atmosphere where the air is denser, 
but not if there are low-lying clouds. 

Notice that we have used our flat-Earth formulas to derive Xdoud, and the conclusions that follow from it. In 
particular, we have expressed the condition for cloud suppression of UASgG rates and (to very good approximation) 
cloud non-suppression of IJAS^S rates analytically without regard to the angles of UAS trajectories. One may 
ask whether Earth’s curvature alters our discussion. Unfortunately, inclusion of curvature leads to transcendental 
equations, rather than to an improved simple analytic expression. However, from numerical studies we can attest 
that curvature does not alter our qualitative conclusions. In fact, when only small horizontal-angle events contribute 
(which holds for most of the cross-section range, cf. Eq. Q), then our quantitative conclusions are accurate, too. 

Finally, we remark that there is no analogue of ^cioud for cloud suppression of HAS rates. The HAS constraints are 
different from the UAS, and HAS trajectories are not restricted to small horizontal-angles. 


V. UAS WITH CLOUDS AND EARTH-CURVATURE 

The distance between interaction point and shower extinction for HAS events is sufficiently small that Earth- 
curvature can be neglected. However, we have seen that Earth-curvature cannot be neglected for UAS events, since 
the tau decay path at high energy provides a length large on the scale of the atmospheric height Zthin- 

We consider again the two UAS possibilities, viewed from space and viewed from ground. The results of sections HV Al 
and llV GI for UAS rates with clouds are extended to include also Earth’s curvature by priming appropriate variables. 

The algorithm for priming was presented in According to the algorithm, the Earth’s curvature is added to 

our prior “flat-Earth” calculations by simply replacing the unprimed variables Zdk, On, 2 :_b(UAS), Z£;(UAS) with primed 
variables z'jj., 0)^, ^^(UAS), z^(UAS), where z'ij.(zdk, On) and 0n{zdk, On) are given in Eqs. (l^ - lt^ . and 2 :^(UAS) and 
0 ^(UAS) are given by priming z_b(UAS), Z£;(UAS), Zdk, S'n in Eqs. and lITTIl to get Eqs. lIHKll and lIHHIl . The 
parameters Zthim -Aground, ^cioud are never primed, for they define layers concentric with the spherical Earth. In what 
follows, we apply this algorithm explicitly to the ground-based and space-based UAS rates. 


A. UAS Viewed from Ground, with Clouds and Curvature 

Priming the appropriate variables of the “too-thin” constraint in Eq. GD becomes 

0 ^ ^ -^^cloud • 

The Zmin constraint Eq. Ga) becomes 

[min{ 2 ;;^(UAS), Zdoud} - Zb(UAS)] tand'n > /min ■ (95) 

This concludes the inclusion of curvature in the integration boundaries of the UASgG rate with clouds. 

To summarize the formulas giving the UASgG events, the general rate equation is Eq. 0 , with inputs from 
Eqs. 0, ®, ®-Ga and (EHll . the “too-thin” constraint is given by (IMIl . and the /min constraint by Eq. TO . with 
z^(UAS) and ^^(UAS) given in Eqs. lIHIill and (IH^ . respectively. 


B. UAS Viewed from Space, with Clouds and Curvature 

Priming appropriately, the “too-thin” constraint in Eq. P?|l becomes 

0 ^ -^dk '^thin ; 


( 96 ) 
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and the Zmin constraint in Eq. becomes 

[min{z^(UAS), Zthin} - niax{ 2 ;B(UAS), Zdoud}] tan6>^ > (97) 

This concludes the inclusion of curvature in the integration boundaries of the UAS^S rate with clouds. 

To summarize the formulas giving the UAS0S events, the general rate equation is Eq. 0, with inputs from 
Eqs. 0, 0, 0-(II21), and P|l . The “too-thin” constraint is given by Eq. and the Zmin constraint by Eq. JiTzli . 
with Zg(UAS) and 2 :^(UAS) given in Eqs. (IHHll and lUTHll . respectively. 


VI. RESULTS 

In this section, we present the results of our semi-analytical approach. For ground- and space-based detectors, we 
show the dependence of the acceptance for HAS and UAS events on neutrino energy, threshold energy, shower length, 
and shower column density, as a function of the neutrino-nucleon cross-section. For incident neutrino energies, we 
choose = 10^° and 10^^ for illustration, and demonstrate that the ratio of HAS-to-UAS events resulting from 
these energies would be of great help in determining the neutrino-nucleon cross-section at these very high energies. 
For the UAS sample, we compute the acceptance for taus emerging over land from pure rock, and separately for taus 
emerging over the ocean from a water layer overlaying a rock layer; we take the water layer to have a uniform depth 
of 3.5 km. We also consider the deleterious effects of low or high cloud layers in the atmosphere, as viewed from space 
and from the ground. 

In al l fig ures, we take the FOV and solid angle entering the acceptance calculations to be that of the EUSO design 
report |35l |. This FOV area, entering Eq. H23II . is tt x (400/'\/3)^ km^. The solid angle is 2 7r for either the HAS or the 
UAS events. The product of area and solid angle is then, very nearly 10® km^ sr. The OWL proposal (two (or 
more) free-flying satellites) has a larger FOV, and stereo eyes. 

Before proceeding with a comparison of the various acceptance curves, it is worthwhile to reflect on what kind of 
event rates might arise in very-large EUSO/OWL-scale neutrino experiments. The event rate is obtained by simply 
multiplying the acceptance by Nature’s cosmic neutrino flux per appropriate flavor. As discussed in earlier sections, 
the appropriate flavor for HAS is i^e (and Df.), since and interactions “lose” 80% of their energy to the escaping 
charged muon or tau. For UAS, the appropriate flavor is Vr, since among the charged leptons only the tau has a 
radiation length long enough to allow a significant fraction of taus to escape from the Earth. For the UAS case, we 
weight the Vr flux by tau branching fractions and r —> shower energy-transfers (2/3 for hadronic showers, 1/3 for 
electronic showers, and zero for the muonic mode). 

The cosmic neutrino flux is a matter for pure speculation at present. A collection of theoretical fluxes is shown in 
Ref. 1^ . We will choose as our benchmark a neutrino flux which is ten times the integrated flux of cosmic-rays at 
Agzk, just below the GZK suppression. Our benchmark (BM) value is 

dJ^BM _ 20 X Eczk) _ 1 

dA dil, dt dA dfl dt km^ sr yr 

The factor of ten is included to give a simple number for the benchmark flux. 

A popular alternative benchmark neutrino flux is that of Waxman and Bahcall 
the high-energy neutrino flux to the observed high-energy cosmic-ray flux. They 

dEwB _ 6 X IQ-^ (10^°eV/U^) 

dA dfl dt km^ sr yr 

Subsequent discussion has shown that their arguments, while sensible, are not compelling. Predictions of the cosmo- 
genic neutrino flux ^3 j resulting from charged-pion production in the GZK process and subsequent pion decay, gives 
fluxes of order of the WB benchmark. Proposed sources of a more exotic nature give larger fluxes. 


(98) 


, who offered arguments relating 
obtained 


A simple estimate of the instantaneous EUSO acceptance for HAS cosmic-ray events is readily obtained by multiplying this Ax 2 tt value 
by 4 to account for the mean projection of the FOV normal to the source. The result is a naive HAS acceptance of ~ 5 X 10® km'^ sr 
for cosmic-rays. For neutrinos, the detection efficiency is less than unity by the factor ~ 2hp{0)a'^^. The factor of two arises because 
the mean path length in the atmosphere of a neutrino is twice the vertical value. Put another way, the increased interaction probability 
for oblique trajectories compensates the ^ coming from projecting the FOV normal to the mean neutrino direction (cosines cancel). 
These simple HAS acceptances assume 100% detection efficiencies. Incidentally, the discriminator between cosmic-ray initiated HAS 
and neutrino initiated HAS is the depth of origin of the shower in the atmosphere. 










[cm^] 



[cm^] 


FIG. 7: Acceptances for space-based (or ground-based) detectors in the absence of clouds. Values of dmin and dmax are fixed 
at 400 and 1200 g/cm^, respectively. The curves correspond to HAS (dotted line), which are independent of Ei, except through 
^’vn{Ev)', and UAS over ocean with = 10^^ eV (thick solid line), ocean with Ei, = 10^° eV (thin solid line), land with 
Ev = 10^^ eV (thick dashed line), and land with Ei, = 10^° eV (thin dashed line). Panels are for (a) StJ) = 10^® eV and 
^min = 10 km; (b) = 10^® eV and Imin = 5 km; (c) E^^ = 5 x 10^® eV and /min = 10 km; (d) Sth = 5 x 10^® eV and /min = 5 

km. For reference, a popular QCD-extrapolation of the neutrino-nucleon cross-section [13 gives 0.54 and 1.2 times 10“®^ cm® 
at Ev = 10®® and 10®^ eV, respectively; the known CC cross-section is 2 x 10“®“^ cm® at an equivalent fixed-target energy of 
5 X 10^® eV, the highest energy for which measurement has been made (at HERA). 


In reality, only Nature knows the value of the real flux. It could be larger than these benchmarks, it could be 
smaller, or it could even be zero. For our benchmark flux, an acceptance of one km^-sr is required to yield one 
event per year. For Nature’s flux, the event rate is ^ / d^diuit ^ ’'^here FI* is the minimum neutrino energy 

producing observable events in the detector with efficiency of order unity. We use units of (km^-sr) when we plot 
acceptances. The fluxes and rates discussed here give a real physical meaning to these acceptance units. 

In Fig.[7|are plotted UAS (solid and dashed) and HAS (dotted) acceptances in our standard units of (km^-sr), versus 
fixed values of cr^^, for the ideal case of a cloudless sky. Five separate dependences are illustrated in this figure: UAS 
vs HAS; El, = 10^® eV vs. 10^^ eV; over ocean vs. over land; shower threshold energy = 10^® eV vs. 5 x 10^® eV; 
and minimum shower length /min = 10 km vs. 5 km. Shower-evolution parameters are set to dmin = 400g/cm® 
and dmax = 1200g/cm®. A sixth possible dependence is whether the shower is viewed from above by a space-based 
observatory, or from the below by a ground-based observatory. Within the approximations of this paper, there is no 
difference between the acceptances for ground-based and space-based detectors in the cloudless case. However, there 
are significant up-down differences when the sky includes clouds. 

The HAS acceptances depend on neutrino energy only via a^^{E,y), and rise linearly with Plotted against 

fixed then, the straight-line HAS curves (dotted) are universal curves valid for any E^, exceeding the trigger 
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threshold The UAS acceptances have a complicated dependence on it arises from the energy dependences 
of V propagation in the Earth, tan propagation in the Earth, and path-length of the tan in the atmosphere before it 
decays, the latter also affecting the visible shower characteristics. In each panel, we show UAS acceptances for two 
different incident neutrino energies, 10^^ eV (thick lines) and 10^° eV (thin lines). The solid lines show the UAS 
acceptances for trajectories emerging from the ocean, and dashed lines show the UAS acceptances for trajectories 
emerging from land, having traveled through only rock. 

Two different shower threshold energies are shown, = 10^® eV in the upper two panels and = 5 x 10^® eV 
in the lower two panels. The EUSO experiment is working to lower its threshold trigger from = 5 x 10^® eV 
to 10^® eV, in order to better overlap events from the Auger experiment (the Auger threshold is ~ 10^® eV). Also 
explored in the different panels of Fig. [71 is the dependence of the acceptance on the minimum shower-length required 
for experimental identification. In the two left panels we have taken = 10 km, while in the two right panels we 
took /niin equal to half of that, 5 km. 

Several trends are evident in Fig.jTj We can clearly see that the UAS acceptance (and so also the rate) is typically 
an order of magnitude larger when neutrinos traverse a layer of ocean water, compared to a trajectory where they 
only cross rock. Thus, the UAS event rate is enhanced over the ocean relative to over land The value of this 
enhancement depends on the shower threshold-energy E®() of the detector (upper versus lower panels) and on the 
neutrino-nucleon cross-section (the abscissa) in a non-trivial way. One sees general trends that (i) the larger the 
cross-section and threshold energy are, the larger is the relative enhancement; (ii) the lower the threshold energy is, 
the closer are the acceptances for different initial neutrino energies (thick vs. thin lines); and for high E®J) approaching 
Ej^, there is a significant suppression of events over land, and over water for larger cross-sections. Lower threshold 
energies are of course also advantageous in that they necessarily imply larger total event rates. 

The sensitivity to E^(j, is partly due to the various energy transfers from the tau to the shower in the different tau- 
decay modes. We have remarked that for the hadronic/electronic/muonic decay modes, |/^/0 of the tau energy goes 
into the shower. This means that a neutrino with an incident energy of 10^® eV characteristically produces a tau with 
energy 0.8 x 10®® eV, which then produces hadronic/electronic showers with mean energies at most (after allowing 
for the tau’s dE/dx in the Earth) 5.3/2.7 x 10^® eV. Clearly, the electronic mode is below the E/j^ = 5 x 10^® eV 
threshold, and the hadronic mode is barely above. Both modes are above the E// = 10^® eV threshold. 

We obtain benchmark event rates by multiplying our calculated acceptances with the benchmark integrated flux, 
Eq. 10811 . of one neutrino per (km® sryr). The result is a signal exceeding an event per year for an acceptance exceeding 
a (km®-sr). Thus we see that the benchmark flux gives a HAS rate exceeding 1/yr if exceeds 10“®® cm®; and 
an UAS rate exceeding 1/yr over water for the whole cross-section range with E^J) = 10^® eV, and over land if 
^ 10“®® cm®. When E^J) is raised to 5 x 10®® eV, however, the UAS signal over land is seriously compromised, 
while UAS rates over the ocean are little changed. HAS rates are unchanged, as long as E// exceeds Ey. 

It is interesting that in the UAS case over the ocean, the acceptances at E^, = 10®® and 10®® eV as a function 
of are seen to cross. For lower values of the cross-section, the acceptance is larger when the initial neutrino 
energy is smaller, unlike what might naively be expected. This is due to the combined effect of larger nadir angles 9n 
contributing at lower neutrino energy, the nature of energy losses in water vs. rock, and the trigger constraints imposed 
on the showers. For larger cross-sections, only small-angle Earth-skimming neutrinos contribute, so the propagation 
of the tau lepton happens mostly in water; complications are mainly absent and hence larger initial neutrino energies 
give larger acceptances, in agreement with intuition. 

We call attention to the fact that for UAS over both ocean and land, there is a maximum in the UAS acceptance 
at cross-section values ~ (1 — 2) x 10“®® cm® and ~ (0.3 — 0.5) x 10“®® cm®, respectively. For cross-sections 
similar or smaller than those at the maximum, the acceptance for UAS is larger than that for HAS; conversely, for 
cross-sections above those at the maximum, HAS events will dominate UAS events. The cross-section value at the 
maximum lies just below the extrapolation of the Standard Model cross-section, which for the two initial neutrino 
energies considered, 10®® eV and 10®® eV, is 0.54 x 10“®®cm® and 1.2 x 10“®®cm®, respectively. If this extrapolation 
is valid, then one would expect comparable acceptances (and event rates) for UAS over water and for HAS; the 
acceptance for UAS over land is down from these by an order of magnitude. If the true cross-section exceeds the 
extrapolation, then HAS events will dominate UAS events; if the true cross-section is suppressed compared to the 
extrapolation, then UAS events will dominate HAS events. Importantly, the very different dependences on the cross- 
section of the HAS (linear) and UAS acceptances offers a practical method to measure One has simply to exploit 
the ratio of UAS-to-HAS event rates. 

Furthermore, the shape of the UAS acceptance with respect to establishes the “can’t lose theorem” P|, which 
states that although a large cross-section is desirable to enhance the HAS rate, a smaller cross-section still provides 
a robust event sample due to the contribution of UAS. The latter is especially true over ocean. 

Finally, from the comparison of left (Zmin = 10 km) and right (/min = 5 km) panels, one infers the sensitivity of 
acceptance to the experimental trigger for visible shower length. We see that reducing the minimum shower length 
by a factor of two here, increases the acceptance by roughly a factor of three for HAS, and slightly less for UAS. 
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FIG. 8: Acceptances in the presence of a cloud layer at «cioud = 2 km; with Zmin fixed at 5 km, dmin at 400g/cm^, and dmax 
at 1200g/cm^. The curves correspond to HAS (dotted line), which are independent of Ei, except through and UAS 

over ocean with E^ = 10^^ eV (thick solid line), ocean with Ev = 10^° eV (thin solid line), land with E^ = 10^^ eV (thick 
dashed line), and land with Ei, = 10^° eV (thin dashed line). Panels are for (a) ground-based detectors with iJtJ) = 10^^ eV; 
(b) spaced-based detectors with = 10^® eV; (c) ground-based detectors with = 5 x 10^® eV; (d) spaced-based detectors 
with E^^ = 5 X 10^^ eV; 


So for a cloudless sky, not too much is lost by choosing longer showers for event reconstruction. This is fortunate, 
for, as remarked early in Section EH the signal/noise and angular reconstruction are greater for longer showers. We 
forewarn that the sensitivity to the ^min trigger will become extreme when we consider a sky with clouds, which we 
address next. 

In Fig. E are shown the acceptances in the presence of a cumulus cloud layer at 2 km, again as a function of cr^^, 
and again with dmin = 400g/cm^ and dmax = 1200g/cm^. We model the cloud layer as infinitely thin with altitude 
-Zcioud, but with an infinite optical depth so that showers are completely hidden on the far side of the cloud layer. 
Details of this modeling were given in Sections CYI and El We call low-lying cloud layers “cumulus”, and high-lying 
layers “cirrus”, for obvious reasons. 

For a sky with clouds, we show acceptances for ^min set to 5 km. We do not show the case with ^min = 10 km, 
because with low-lying clouds, the UAS acceptances for ground-based detectors are essentially zero with ^min = 10 km 
(and dmin = 400g/cm^ and dmax = 1200g/cm^). On the other hand, the space-based rates are virtually unaffected 
by the low clouds. For the UAS as seen from the ground, there simply is not enough space below the cloud layer for 
the Earth-skimming tau to decay and for the subsequent shower to develop (see Eqs. Gi and (IH3)- The smaller 
^min = 5 km that we do show allows enough UAS events to develop into an observable shower below the cloud layer 
to establish a meaningful acceptance for ground-based detectors. Recall, however, that the UAS rate suppression 
due to clouds depends sensitively on the value of (dmax — dmin)/^min- If the value of this is larger than dvert/h, the 
suppression is aggravated for space-based UAS and alleviated for ground-based UAS, and vice versa (as discussed 













FIG. 9: Dependence of acceptance on cloud altitudes for space-based fluorescence detectors. Fixed values are Zmin = 5 km, and 
threshold energies of = 10*^® eV in the left panel and 5 x 10^® eV in the right panel. All curves representing UAS assume 
trajectories over water and an initial neutrino energy of 10®® eV; curves for HAS are valid for any energy exceeding i?th- Solid 
lines show HAS, while dashed and dotted lines show UAS with and without Earth-curvature effects, respectively. Thick lines 
are for a cloud layer at 4 km and thin lines for a cloud layer at 12 km. 


below Eqs. m and m, and more extensively in Section ITV ED . 

Since we show one value of Zmin, not two, in Fig. |H1 the number of panels is half that in Fig. [3 On the other 
hand, the symmetry between upward-looking ground-based detectors, and downward-looking space-based detectors 
is broken by the cloud layer, so we must now show separate panels for the space-based and ground-based detectors. 
This brings the number of panels back to four. The HAS and UAS curves are represented in the same way in Fig.|3 
as for Fig. [3 The thick and thin lines bear the same meanings for the initial neutrino energies. The left panels 
show acceptances for ground-based detectors, whereas the right panels show those of space-based detectors. As with 
the previous figure, the upper panels show results for a threshold energy of = 10^® eV, and the lower panels for 
Ftt = 5 X 10^® eV. 

We see from this figure that qualitative features learned for the cloudless case apply also in this cloudy case. One 
difference is that the UAS acceptances over land in panel (c) are smaller than the HAS acceptance over the entire 
range of cr®^. One may judge the effect of clouds by comparing Fig.|51against the Zmin = 5 km panels (b) and (d) of the 
clear-sky Fig. [3 Quantitatively, the ground-based acceptances (left panels in Fig. noclouds) are quite reduced by the 
low-lying clouds, whereas the space-based acceptances (right panels) are not, as one would expect. The suppression of 
the ground-based acceptance is most severe for small cross-sections, for which the tau leptons emerge more vertically 
and disappear into the clouds before their eventual shower occurs and develops. Ground-based UAS acceptances are 
reduced by up to an order of magnitude over water, and even more over land. Ground-based HAS acceptances, still 
linear in cr®^, are reduced by an order of magnitude. For space-based detectors, the UAS acceptance is reduced 
little by clouds at 2 km. Larger neutrino cross-sections lead to more tangential tau-showers which may hide below a 
low-lying cloud layer. We see that UAS reductions are a factor of 2 for the larger cross-sections shown, and less for 
the smaller values of cross section. 

The dramatic reduction of ground-based acceptances by low-lying cumulus clouds begs the question, “what are the 
effects of higher-altitude clouds on space-based detectors?” In Fig. |3 we continue the study of the dependence of 
space-based acceptances on cloud altitude. We also examine the suppressing effect of the Earth’s curvature. 

Fixed values in Fig. |3 are Zmin = 5 km, and threshold energies of E’jJ) = 10^® eV in the left panel and 5 x 10®® eV 
in the right panel. All curves representing UAS assume trajectories over water and an initial neutrino energy of 
10®® eV; curves for HAS are valid for any energy exceeding Acceptances for two cloud altitudes, Zcioud = 4 km 
(thick curves) and 12 km (thin curves), are shown for HAS (solid curves) and UAS (dashed curves). Results are 
to be compared with the thin solid (UAS) and thin dotted (HAS) lines in panels (b) and (d) of Figs. Hand|3 We 
infer from comparing these three figures that the effect on a space-based detector of higher cumulus clouds, and even 








FIG. 10: Dependences on dmin, dmax, and Zmin without clouds of acceptances for space-based fluorescence detectors. Thresh¬ 
old energies are = 10^® eV in the left panel, and 5 x 10*^® eV in the right panel. All curves representing UAS assume 
trajectories over water, and an initial neutrino energy of 10®® eV; curves for HAS are valid for any energy exceeding i?th- 
Thick lines represent HAS, and thin lines UAS. Solid lines correspond to (dmin, dmax, ^min) = (Og/cm®, 10® g/cm®, 0 km), 
dashed lines to (300g/cm®, 1500g/cm®, 5 km), dotted lines to (300 g/cm®, 1500 g/cm®, 10 km), and dash-dotted lines to 
(400g/cm®, 1200g/cm®, 10 km). 


higher cirrus clouds, is more dramatic for downgoing HAS than for upcoming UAS. The HAS acceptance is reduced 
by factors of ~ 1.5, 3, and 10 when the cloud layer lies at Zcioud = 2 km, 4 km, and 12 km, respectively. In contrast, 
the UAS acceptance is reduced by factors of ~ 1.5, 2 and 3 when the cloud layer lies at Zcioud = 2 km, 4 km and 
12 km, respectively. Since cloud layers are common, they will compromise the acceptance of space-based detectors. 

Also shown in both panels of Fig. 0 are the UAS acceptances (dotted lines) for a flat Earth. One sees that correct 
inclusion of the Earth’s curvature lowers the acceptance, since it puts the tau decay and the subsequent onset of 
shower evolution into the thinner air of higher altitudes. Curvature does little harm for smaller cross-sections, but 
reduces the acceptance for ^ 0.5 x 10“^^ cm^. Coincidentally, 0.5 x 10“^^ cm^ is the popular value for the QCD- 
extrapolated cross-section. The reduction of acceptance for for larger cross-sections is understandable, because the 
for larger the taus emerge from the Earth more horizontally, and hence travel more lateral distance before they 
decay. The Earth “falls away” from the taus as (lateral displacement)^/2 Beyond ~ 10“^^ cm^, the reduction 
factor is about 2.5 for cloud layers at either 4 or 12 km (and quite different for E^, near E^^J. 

Since curvature raises the altitude of the tau-shower (and rotates it toward the vertical by 6 * ^ (lateral 
displacement)/i? 0 ), the net effect is to remove the bottom layer of atmosphere from the UAS shower development. 
Clouds remove the bottom layer from view for UASCiS. Thus, one expects the reduction in acceptance due to Earth’s 
curvature to be largest in the cloudless case. We have checked numerically that for the parameters of Fig. O a 
reduction factor of ~ 3.5 is obtained for the cloudless case. 

However, it is dangerous to generalize that Earth’s curvature causes event suppression. One sees in the left panel 
of Fig. El that for high clouds, curvature effects may even increase the event rates. This is because more events that 
would not have been visible above the cloud altitude are now “lifted” to visibility, compared to the number of events 
that would have been visible but are not “lifted” to invisibility. Curvature does not increase rates in the right panel, 
which points again to the dangers of generalization. 

In Fig.^1 we show dependences of the acceptances on the parameters Zmin and dmin describing the shower triggers, 
and the parameter dmax characterizing shower extinction. In this figure we assume a cloudless sky. Threshold energies 
are = 10^® eV in the left panel and 5 x 10^® eV in the right panel. Thick curves present HAS acceptances and thin 
curves present UAS acceptances. All curves representing UAS assume trajectories over water, and an initial neutrino 
energy of 10®® eV; curves for HAS are valid for any energy exceeding We depict four different sets of shower 
parameters: (dmin, dmax, ^min) = (0g/cm^, 10® g/cm^, 0 km) shown in solid curves, (300g/cm^, 1500g/cm^,5 km) in 
dashed curves, (300g/cm^, 1500g/cm^, 10 km) in dotted curves, and (400g/cm^, 1200g/cm^, 10 km) in dash-dotted 
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curves. The value dmax = 10® g/cm^ is to be interpreted as an effectively infinitely long shower persistence, i.e. an 
illustration of shower development without extinction. The acceptances in the solid curves correspond to the most 
liberal shower-trigger requirements, with basically all showers declared observable. Those in he dashed, dotted,and 
dashed-dotted curves correspond to a realistic set of choices for the shower-development parameters dmin, dmax and 
for the shower length /min- 

One sees that even in the cloudless case shown here, the dependences on the shower parameters is considerable. 
The UAS acceptance gets reduced by ~ 2 when Zmin is increased from 5 km to 10 km and the visible column density 
(dmax — dmin) IS reduced from 1200 to 800 g/cm^. Changes in the HAS acceptance are a bit more dramatic. As the 
shower triggers are tightened according to our examples, the HAS acceptance falls by a factor ~ 6. 

The acceptances in Fig.^lcan also be compared to the thin solid (UAS) and dotted (HAS) curves in panels (b) and 
(d) of the no-cloud Fig.[7| where (dmin, dmax, ^min) = (400g/cm^, 1200 g/cm^, 5 km). As can be seen from this figures, 
the effect of reducing the constraint on the minimum shower length from 10 km to 5 km increases the acceptance by 
roughly a factor of ^ 3 for HAS, and ^ 2 for UAS. 

With clouds, the sensitivity to shower-development parameters is more acute. The dotted case corresponds to 
a good trigger sensitivity (discussed in Section HV Fll of (dmax — dmin)/^min = 120gcm“^/km, and so to a critical 
cloud-altitude of = 0.56 km (refer to Fig.lHJl. The dashed-dotted case corresponds to a less good 80gcm“^/km 
sensitivity, and a less pleasing (for ground-based detectors) critical cloud-altitude = 3.8 km. Thus one expects 
the dashed-dotted case to be quite sensitive to cloud layers, and the other three cases to be relatively insensitive to 
cloud layers. 


VII. DISCUSSION 

In this paper we have presented a mostly analytic calculation of the acceptances of space-based and ground-based 
fluorescence detectors of air-showers at extreme-energies. Included in the calculation are the dependences of the 
acceptances on initial neutrino energy, trigger-threshold for the shower energy, composition of Earth (surface rock 
or ocean water), and several shower parameters (the minimum and maximum column densities for shower visibility, 
and the tangent length of the shower). Also included in the calculation are suppression of the acceptances by cloud 
layers of arbitrary altitude, and in the UAS case, by the Earth’s curvature. And most importantly, included in the 
calculations are the dependences on the unknown neutrino cross section. The dependence is trivial and linear for 
HAS, but nontrivial and nonlinear for UAS. 

The merits of the analytic construction are two-fold: it offers an intuitive understanding of each ingredient entering 
the total calculation; and it allows one to easily re-compute when parameters governing the shower or the atmosphere 
are varied. While a Monte Carlo approach may be simpler to implement, it sacrifices some insight and efficiency. 

The differing dependences of HAS and UAS on enable two very positive conclusions: (1) the “no-lose theorem” 
is valid, namely, that acceptances are robust for the combined HAS plus UAS signal regardless of the cross-section 
value; and (2) an inference of the cross-section at 10^° eV is possible if HAS and UAS are both measured. 

Our formulas are valid for the energy range ^ 10^® —10^^ eV. The lower limit is necessary to validate our assumption 
that the r decay length is much larger than its radiation length. Below 10®® eV, the t’s boost factor is insufficient 
to provide a t lifetime in compliance with this assumption. The upper limit arises from the fact that above 10^® eV, 
the weak charged-current losses of the t’s, not included in our calculation, exceed their electromagnetic losses, which 
we have included. Although the energy range of validity is limited, it covers the range of interest for extreme-energy 
cosmic neutrino studies. 


A. A tale of two media (ocean and land), and two lengths (Xv and Ar) 


One may sensibly ask, “what difference does it make whether the UAS interaction takes place in water or in 
rock?” After all, one has only to look at a different value of L or to compensate for any change in target density. 
Nevertheless, we have shown that the UAS rate over ocean is an order of magnitude larger than over land. Let us 
explore this a bit. 

The rate for producing upgoing taus which exit the Earth’s surface has its peak near the chord length Apeak ~ X^+Xt- 
Increasing the solid angle optimizes the rate. This is accomplished by making Apeak as large as possible. Equivalently 
then, we ask that Xi, and be as large as possible. Explicitly, we seek to maximize 


A 


peak ‘ 



Xu + A- 


Pearth 


Na /?. 


( 100 ) 
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Apparently, Nature may optimize Apeak in one of three ways: either reducing the cross-section, reducing the tau’s 
radiation loss, or reducing the density. Thus, the oceans, less dense than surface rock, will serve to increase the 
emerging tau rate by the factor psi/pw = 2.65 (for fixed and /?,-)■ For generation of the visible UAS, there is a 
further enhancement. Trajectories through water exit the Earth nearly horizontally. Horizontally-emergent taus have 
more time and more atmosphere in which to decay and evolve into a visible shower. This further increases the UAS 
acceptance. And there is a third effect: Pt in water is nearly half of Pt in rock, leading to a larger Ar in Eq. (trnnii . 
So several effects conspire to enhance the UAS acceptance over water vs. over land. In our numerical work presented 
in Figs.[7|andO we found the enhancement to be typically a factor of 10. This enhancement over water is fortunate, 
in that 70% of the Earth’s surface is ocean, and ocean nights are not polluted with man-made lights. 

Note that 


Ar _ / Pig \ ^ 1 0-06 for rock 

Ay \Pr{E)) I 0.10 for water 


( 101 ) 


Furthermore, the energy dependence of Pt{E) oc is mild, changing P^ by just 1.6 per decade. Thus, we have that 
Ay > Ar for ^ 10“^° cm^, i.e. for the entire range of cross-sections which we consider. So the story of the UAS 
enhancement over water is really the story of the neutrino’s longer MFP in water than in rock, followed by the tau’s 
greater probability to decay and shower following its nearly tangential emergence from water. 

For very we ak cross-se ctions, the MFP becomes large compared to the size of the critical chord in the ocean 
Aw-sr = 422 •^/^((TSAlcm km. Then the neutrino is forced to spend part of its trajectory in rock, mitigating the 
difference between UAS over water and over land. This phenomenon is included in our calculations. We have seen 
that the neutrino MFP exceeds the critical chord for < 4 x 10“^^cm^. Far below this cross-section value, the 
ocean portion of the neutrino’s chord is too small to affect the UAS rate. In Fig. [7| one can see that the two UAS 
acceptances, over water and over land, tend toward each other below ^ 10 “^^ cm^, and to a ten-fold enhancement 
for water above ~ 10 “^^ cm^. 


B. Remark on Observations over Elevated Land 

We will present results for events over land and ocean, both taken to have zero, i.e. “sea level”, elevation. However, 
much of the Earth’s land surface is at higher elevation (fortunately, for land animals). Furthermore, some ground- 
based observatories are sited at high elevations to reduce various backgrounds. The Auger observatory in Argentina, 
for example, is at 1400 meters, above ~ 15% of the atmosphere. The HiRes siting in Utah is at a similar elevation. 
For HAS viewed from these sites (in cloudless skies), the acceptance is the same as that for space-based viewing with 
a cloud layer at 1.4 km. For UAS viewed from elevated sites, one cannot proceed by simple analogy. 

However, inclusion of elevation into our formalism, for ground-based or space-based, HAS or UAS, with clouds and 
without, is simple. One merely replaces the sea-level atmospheric density, /o(0), with the ground-level density. For an 
elevation of Zeiev, the replacement value is p(0)- This leads to the further replacement dvert ^ dvert- 

With these substitutions, all previous formulas may be used, with altitudes of clouds and other “z-parameters” now 
understood to be with respect to ground-level, not sea-level. Strictly speaking, zthin should be replaced by zthin —-^eiev 
This would allow accurate comparisons of acceptances from elevation with those from sea-level. In practice, keeping or 
not keeping the additional kilometer or so of thin atmosphere makes little difference. And for comparing ground-based 
and space-based detectors at a common elevation, the difference is irrelevant. 

In summary, ground elevation reduces the amount of atmospheric volume available, sometimes substantially. This in 
turn reduces the target mass for HAS, the decay volume for UAS, and the grammage available for shower development. 
Consequently, elevation disadvantages observations over land, compared to space-based observations over the zero- 
elevation ocean. The disparity between the two becomes more acute in the presence of clouds. 


C. Comparisons with prior work 

Among the aims of this paper is the detailed extension of the idea introduced in Ref. Q , that at 10^° eV can be 
inferred from a measurement of the UAS/HAS ratio. So let us first compare our calculation with the more qualitative 
one in Ref. P|. In Ref. Q, several energy-dependences were purposely frozen, for simplicity. The tau energy loss 
was set constant Pt{E) = Pig (i.e., a = 0), and the energy for the produced tau lepton was assumed to be that of 
the incoming neutrino (i.e., (y) = 0). The decay length for the taus exiting the Earth was fixed to a single value, 
ignoring the dependence on the tau’s initial energy, production point in the Earth, and zenith angle. Furthermore, 
all taus were considered to decay within an atmospheric height of 10 km. These approximations do not affect the 
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main conclusions of P|, but do impact the results quantitatively. For example, they lead to unphysical behavior in 
the UAS rate for high values of the cross-section; for large the incident neutrino must interact with the Earth’s 
surface with probability one, and the UAS acceptance should asymptote to a cr^^-independent constant determined 
solely by tau physics. Also, the case where the trajectory’s chord length in the Earth is smaller than the distance 
for which the tau lepton energy is reduced to the threshold value is incorrectly calculated. The impact of this 
is small, for it only affects very Earth-skimming neutrinos (6*hor < 0.1°). Einally, Ref. only considered neutrinos 
traveling through rock, not water, and did not include suppression effects from realistic shower formation, clouds, or 
curvature of the Earth. Our work considerably improves upon the original work of Ref. Q . 

The prior work in Refs. (and related work in Ref. is semi-analytic, like our own. A main difference 

between them and us is the manner in which the tau energy loss is parametrized. The parameter is taken to be 
constant in Ref. ^3 (where (y) = 0 is also assumed), to depend line arly on the tau energy in Ref. [33, and assigned 
a logarithmic dependence on the energy in Ref. [sj. Also, Refs. |U|43 assume some maximum tau decay distance, 
and do not implement any constraints from the subsequent shower formation. Furthermore, they do not consider 
events over water, or the presence of clouds. Ref. computes the flux of tau leptons exiting the Earth, but doesn’t 
consider the important process of tau decay and shower formation. On the other hand, does include the possibility 
of the tau decaying inside the Earth. However, as we mentioned earlier in our paper, including tau decay within the 
Earth reduces the UAS acceptance by less than a few per cent for the energies of interest here. Also in Ref. a 
Monte Carlo calculation is performed and shown to agree with the semi-analytical approach to very good accuracy. 
Where comparisons are possible, our results agree qualitatively with . Let us note that Ref. offers an improved 
and more detailed evaluation of the effective acceptance for fluorescence detection at the Pierre Auger Observatory. 
Taken into account is the real elevation profile of the surrounding mountains. They find a significant increase in the 
event rate due to the nearby mountains, compared to the semi-analytical, mountain-less calculations of Refs. |33l l42| 
and us. Unfortunately, even with the enhancement, the predicted event rates are ^ 0.5 events/yr for neutrino fluxes 
motivated by the GZK process and topological defect decay models. This small rate points again to one of the major 
benefits of space-based detectors: the much larger FOV. 

In Ref. 1^, a detailed semi-analytical computation of UAS and HAS was performed. This work considered a larger 
range of energies than we do. Hence, it was necessary for Ref. to include the possibility of tau decay inside 
the Earth. This work assumed (3r to be constant, but otherwise the energy-loss of the produced tau was calculated 
accurately. In addition, constraints on the shower formation were included in a semi-analytical approach, in order to 
calculate event rates for an air shower array. The main differences between this work and our calculation is that we 
consider a water layer as well as rock, and we include the possibility of clouds. 

Detailed Monte Carlo simulations of acceptances are presented in Refs. Neutrino scattering inelasticities 

and tau energy losses are accurately included. We have checked that we get very good agreement with the results of 
these papers. Ref. is specific ground-based Auger detector. It includes realistic shower formation and detector 
response, but it does not consider clouds. It also considers only neutrino and tau propagation in rock. While rock is 
the dominant material in the vicinity of Auger, there are trajectories reaching Auger from the West which will travel 
in the Pacific Ocean. Ref. [^. undertaken mainly with EUSO in mind, does study acceptances over both water and 
land. 

In Ref. 1^. the analytic approach is different, and more optimistic rates are obtained. However, there are several 
questionable approximations. How the energy threshold constraint is implemented is obscure. The treatment of the 
Earth’s atmosphere is too simplistic. And although the calculation is meant to be valid for an arbitrary Earth density 
profile, the derived expression for the event number is only valid for a constant Earth density. Our acceptances for 
UAS events do not support the optimistic UAS rates of Ref. ^3- We do however support the results found by the 
more detailed analysis, e.g. that in Ref. [i^. 

In summary, the main advances we present in our study are the inclusion of a new analytic and accurate power-law 
parametrization for the tau energy loss pre-staging UAS events, the analytical implementation of shower constraints, 
cloud boundaries, and Earth curvature, and the consideration of UAS events over the ocean as well as over land. We 
note that the tau energy loss parametrization we implement in the present study was already used in Ref. |3n| for the 
case of taus propagating in rock. Here we present also the parametrization when taus cross a water layer. 


D. Odds and Ends 


In this subsection we offer remarks on issues possibly relevant to this paper. 
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1. Incident neutrino flavor ratios 

One of the main points of this paper is to hone the argument that the CC neutrino cross-section can be inferred 
from a comparison of UAS and HAS rates. The ratio of these rates is the product of a flux ratio times acceptance 
ratio. We have focused on the electron-neutrino as the primary particle for HAS initiation, and the tau-neutrino as 
the primary particle for UAS initiation, for the good reasons given in the text. The relevant flux ratio, therefore, is 
the ratio of the flux to the Vr flux. 

We have calculated acceptances. From these, one may simply form the acceptance ratio. What is not known at 
present is the relevant to i^t flux ratio at ~ 10^° eV. A general theorem for neutrino flavor-mixing states that if the 
atmospheric mixing angle is nearly maximal (it is!), and if the short-baseline angle O 13 is nearly zero (it is!), then 
and i>T equilibrate over cosmic distances. A corollary to the theorem then, is that if cosmic neutrinos originate from 
the complete pion decay chain, TT—ffj, + Vfj,—^e + i>e + then after equilibration the neutrinos will arrive at Earth 
with the democratic flavor ratio of 1:1:1. However, dynamics at the source, or new physics enroute from the source, 
could alter this favorable ratio. Caveat emptorl 


2. Ratio of neutral- and eharged-current eross-sections 

It is an implicit assumption in this work that the ratio of the neutral to charged-current cross-section is small, 
~ 0.44 according to the Standard Model of particle physics. However, it is possible that above 10^® eV and below 
10^° eV a threshold is passed at which the NC interaction becomes strong and the CC interaction does not. Such 
would be the case, for example, in models of low-scale gravity unification. 

Crossing such a hypothetical threshold would change the physics in this paper dramatically. First of all, even 
though the NC interaction typically puts ~ 5 times less energy into the shower than does the Ve CC interaction, with 
a much larger NC cross-section, even at fixed the NC events would dominate the I'e CC events. Secondly, UAS 
acceptances would be reduced because the energy losses of neutrinos passing through the Earth would be larger. 


3. Weakly-interacting non-neutrino primaries 

The range of cross-sections we consider in this work spans the cross sections of weakly-interacting massive particles 
(WIMP), a popular candidate for dark matter. We believe, therefore, that our figures may be useful in assessing the 
qualitative features of acceptances for WIMP detection. However, we caution that there are substantial differences 
between WIMP-initiation of showers and neutrino-initiation. The WIMP carries considerable inertial mass, and so 
transfers less energy to its shower. Also, the UAS generated by a WIMP flux would not proceed through the tau 
production and decay chain. 


VIII. CONCLUSIONS 

We have presented analytic formulas for the acceptances of fluorescence detectors, both space-based and ground- 
based, for neutrino-initiated events, as a function of the unknown extreme-energy neutrino cross-section. For the 
downgoing HAS events, the dependence of acceptance on cross-section is linear, but for upcoming UAS events the 
acceptance is quite complicated. It turns out to be somewhat flat and relatively large, which validates the “can’t lose” 
theorem which says that if the HAS rate is suppressed by a small then the UAS rate compensates to establish 
a robust signal. 

We have studied the dependence of acceptances on the incident neutrino energy, the trigger-energy for the 
shower, shower-development parameters dmin and dmaxi and observable (tangent) shower length Zmin; and on the 
“environmental” conditions of cloud layers for HAS and UAS, and events over ocean versus over land for UAS. UAS 
showers typically originate at a considerable distance {cTr = 4900 (At/IO^® eV) km) from the point on the Earth where 
the parent tau emerged. Therefore, due to the earth’s curvature, they originate at higher altitudes with thinner air. 
Thus, it is necessary to include the Earth’s curvature in the calculation of UAS acceptances. We have done so. We find 
that inclusion of the Earth’s curvature reduces the UAS acceptance by a factor of a few when ^ 0.5 x 10“^^ cm^. 
The meaning of “a few” depends on the various parameters entering the calculation. 

Clearly lower shower-trigger energies are better. This is especially true when clouds are present. We have quantified 
the sensitivity to by comparing two realistic values, 10^® eV and 5 x 10^® eV in the face of incident neutrino 
energies of I0^° eV and 10^^ eV. 
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Cloud layers may severely suppress acceptances. For UAS acceptances, there is a strong dependence on the combi¬ 
nation of shower-trigger parameters (dmax ~ di„in)/^min) especially with clouds present. Maximizing this combination 
to a value of dvert/h = /o(0) = 129gcm“^/km or greater, significantly minimizes the suppression from clouds. For 
(dmax-drain)/^min ^ p(0), there is a critical altitude Zcioud — [(dmax - drain)/^rain) (h/dvert)] below which a cloud 

layer would totally obscure the acceptance of a ground-based detector, but leave the acceptance of a space-based 
detector unaltered. Clouds above the critical altitude would partially obscure UAS events, and therefore suppress the 
acceptances, of both space-based and ground-based detectors. 

Concerning UAS events over water versus over land, we find that acceptances over water are larger, typically by 
an order of magnitude. We have traced this enhancement over water to the increased pathlength in water of both 
neutrinos and taus, and to the increased pathlength in air for tau decay and increased column density in air for shower 
development, when a tau emerges with small horizontal angle from the relatively shallow ocean. We also noted the 
smaller enhancement from the fact that the atmospheric grammage over water integrates from sea-level, whereas 
the grammage over land is often 15% less. It is difficult to imagine a ground-based detector over the ocean, so the 
“water advantage” clearly belongs to the orbiting space-based detectors. Perhaps a ground based detector could be 
positioned near an ocean to realize the “water advantage” for much of its solid angle. 

In the spirit with which we began this study, we are led to two bottom-line conclusions: 

(i) Inference of the neutrino cross-section at and above 10^° eV from the ratio of UAS and HAS events appears 
feasible, assuming that a neutrino flux exists at these energies. 

(ii) Space-based detectors enjoy advantages over ground-based detectors for enhancing the event rate. The advantages 
are a much higher UAS rate over water compared to land, and the obvious advantage that space-based FOV’s greatly 
exceed ground-based FOV’s. 

Our hope is that space-based fluorescence-detection becomes a reality, so that the advantages of point (ii) can be used 
to discover/explore the extreme-energy cosmic neutrino flux. According to point (i), part of the discovery/exploration 
can be the inference of the neutrino cross-section at E^, ~ 10^° eV. 
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APPENDIX A: DENSITY INTEGRATIONS IN THE EARTH 

Two column densities in the Earth are of relevance for UAS probabilities. The first is the path-integral of Pearth for 
the incident neutrino from entrance to interaction in the Earth. This column density controls the neutrino absorption 
probability, and therefore, the neutrino survival probability to the point of interaction Wint- The column density 
dv{0n) is stated in Eq. ® as 

dv{0n) = / die Pearth (w) . (Al) 

For constant density, which applies only for Earth-skimming neutrinos entirely in surface rock (pearth = Psr = 
2.65 g/cm^) or entirely in ocean water (pearth = Pw = 1-Og/cm^), the result is simply = {L — u>int) Pearth ■ 

This constant density result holds in ocean for angles relative to the horizon smaller than 1.90°, and it holds in rock 
for angles relative to the horizon smaller than 22.17°, as we show below. 

The second relevant column density in the Earth is that of the emerging tau in the UAS event sequence. This 
column density dr {On) is the path-integral of Pearth from the interaction point in the earth to the earth’s surface, 

dT{0z)= / dw . (A2) 

Eor constant density, which applies for taus emerging from rock, or for taus emerging from water with horizon-angle 
less than 1.90°, the result is simply dr(dz) = 'Wint Pearth • This column density, when suitably weighted with the tau 
energy-attenuation factor /Jig, controls the tau energy-loss probability, and therefore, the probably for the tau energy 
and lifetime at emergence from the Earth as given in Eqs. 0 and 1I3- 
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Boundary 

Radius 

Critical L 

Critical On 

Critical Ohor 

Earth/atmosphere 

A® = 6371 km 

0 km 

90° 

0° 

water/surface-rock 

rsr = A© — Zw, Zw = 3.5 km 

Lw-sr — 422 ( 3 km 

88.10° 

i-so" (^)^ 

surface-rock/mantle 

Tm = 5900 km 

Lsr-m = 4808 km 

67.83° 

22.17° 

mantle/core 

rc = 3485.7 km 

Lm-c = 10666 km 

33.17° 

56.83° 


TABLE IV: Critical chord lengths and nadir and horizon angles for the indicated Earth bonndaries. 


For our purposes, concentric shells of constant density provide a sufficiently accurate approximation to the Earth’s 
profile, and allow for an analytic evaluation of the path integrals. We take a simple model of this kind for the Earth 
density: 


Pearth(r ) 


Pw = l- 0 g/cm^ for 
Psr = 2.65 g/cTc? for 
Pni = 4.0g/cm^ for 
Pc = 12 . 0 g/cm^ for 


fsr — ■^w ^ r ^ , 

rm <r <rsr, 

Tc < r < Tm , 

0 < r < Tc , 


(A3) 


with the radii of the boundaries listed in Table (©• For UAS events over land, we replace the outermost Zv, = 3.5 km 
of water with surface rock. Thus, there are in this Earth model four (three) concentric density zones for UAS events 
over water (land). 

The chord length L, nadir angle and horizon angle 0hor are related to the sagitta s by the formulas 


L{s) 
cos 9n{s) 
sin On{s) 


and to the boundary radius tb by 


2 'v/ 2 i ?0 s - , 

g 

COs6lhor(s) = 1 - ^ ; 


(A4) 

(A5) 

(A 6 ) 


M^b) 

II 

to 

1 

(A7) 

cos On(rB) 

= si„»Ho,(rB) = , 

(AS) 

sin On(rB) 

r’B 

= COs6»hor(?’B) = ^ . 

(A9) 


In Table El we collect the critical values for L and at the various boundary layers. 


A further useful formula is the path-length w{rB] L) from the Earth’s surface to the boundary of radius tb, for 
fixed L ov On'. 


w{rs; L or On) = —- 



^.2 _ 7?2 


= COS On - 


— Rts, sin 


(AlO) 

(All) 


Consider the calculation of dtot{L). Using the path-lengths defined in Eq. irmii . the result is: 


^tot - 


L , for 0 L -^w—sr 

2 L') Pw “t” ^')] Psr} ; ^w—sr ^ -^sr—m 

2 {^/^(rsr^ L') L') w(Tst'i -^)] Psr ["2 -^)] Pm} ; fol" -^sr—m ^ ^ ^m —c 

2 {rt;(rsr) pw + [r(;(rm) - '^/^(rsr)] Psr + [w{rc) - r(;(rm)] Pm + [f - Pc} , for Lm-c < L < 2R^ 

(A12) 
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where for compactness, we have suppressed the L-dependence of the ui-function in the final line. For events over land, 
Pw must be replaced in Eq. llXnil with psi-- 

Next we consider the calculation of dr{L). The tau radiation length is very short on the scale of r„i, and so the tau 
path is confined to just surface rock and ocean. Over land, then, we have simply (iT(land) = WintPsr- Over water, the 
calculation has two contributions in general, from water and from surface-rock. For L < Lw-sd the tau encounters 
just water, and so c?r(ocean; L < Lw-sr) = WintPw- For L > Lw-sr, the tau encounters rock and then water. However, 
the tau never encounters first water and then rock and then water again, for this requires a tau trajectory exceeding 

1/2 of the critical path-length Lw-sr, which is 211 ^ km, far exceeding the tau radiation length. We summarize 

these results, again making use of Eq. ITmll : 


{ Wint Psr , over land, for all Wint 

Wint Pw , over oceans, for Wint < w{rsr', L) (A-13) 

wiTsr] L) Pw + [wint “ w{rsr] L)] Psr , over oceans, for Wint > w{rsP, L) 

To obtain we may use the simple relation d^, = dtot — dr- Thus, we are finished with calculating column densities. 

The tau opacity is easily obtained in the constant-density concentric-shells approximation. Weighting the segments 
in Eq. ifm with the corresponding values of /3i9, either /3fg = 1.0 x 10 ® cm^/g or = 0.55 x 10 ® cm^/g, we 
have 


Pfg Psr , 
Wint PTg Pw , 


over land, for all Wint 
over oceans, for Wint < 'w{rsi; L) 
w{rsi; L) /3^g pw -I- [wint - w{rsi;L)\ psi , over oceans, for Wint > w{rsi; L) 


(A14) 


We also need an explicit formula for wth(i), defined implicitly in Eqs. @ and m- For notational brevity, let us 
recall the notation in Eq. ca: 


I{wth) = 


10i®eV 


th 


10i®eV 

0.8 L,, 


(A15) 


Then, a calculation similar to the ones above leads to 


{ 7 ^— X{wth ), over land, for all L , 

Pl9 psr 

I{wth), over oceans, for L < Lw-sr or I{wth) < PTg Pw w(?’sr; L) , 

w{rsi; L) (^1 - ^ -I- T(wth), over oceans, for L > Lw-sr and l{wth) > Pw w(?’sr; L). 

(A16) 

Finally we note that, in all the above formulas for events over oceans, the correct result over land may be found by 
simply setting pw equal to psr and Pig equal to /3fg. 
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